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Abstract. We study F-convergence of graph based Ginzburg-Landau functionals, both the hmit for zero 
diffusive interface parameter e — >■ and the limit for infinite nodes in the graph m — >■ oo. For general graphs 
we prove that in the limit e — >■ the graph cut objective function is recovered. We show that the continuum 
limit of this objective function on 4-regular graphs is related to the total variation seminorm and compare 
it with the limit of the discretized Ginzburg-Landau functional. For both functionals we also study the 
simultaneous limit e — )■ and m — >■ oo, by expressing e as a power of m and taking m — ^ oo. Finally we 
investigate the continuum limit for a nonlocal means type functional on a completely connected graph. 



AMS Subject Classifications: 35R02, 35Q56 

1 Introduction 

1.1 The continuum Ginzburg-Landau functional 

In this paper we study an adaptation of the classical real Ginzburg-Landau (also called Allen-Cahn) 
functional to graphs. The Ginzburg-Landau functional is the object to be minimizec^ in a well known 
phase field model for phase separation in materials science, e.g. |451 H6] and is given by 



where u £ W^''^{Q) is the phase field describing the different phases the material can be in and is a 
double well potential with two minima, e.g. W{s) = — 1)^. is a bounded domain in M^. 

Recently [HI this functional has been adapted to weighted graphs in an application to machine 
learning and data clustering: An image is interpreted as a weighted graph, with the vertices corre- 
sponding to the pixels and the weights based on the similarities between the pixels' neighborhoods. 

^Note that to avoid trivial minimizers an additional constraint needs to be added. In materials science it is common 
to add a mass constraint of the form f^u — M for a fixed M > 0. In image analysis applications one often adds a fidelity 
term of the form X\\u — /||^2(q) to the functional F^^, where A > is a parameter and / £ L^(Jl) is given data, often a 
noisy image which needs to be cleaned up, [52) . 
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The phase separating nature of the Ginzburg-Landau functional then drives separation of the different 
features in the image. 

The continuum functional F^^ has been extensively used and studied, but a theoretical under- 
standing of its equivalent on graphs is lacking. In this paper we use T-convergence [281 US] to study 
the asymptotic behavior of minimizers of the graph based Ginzburg-Landau functional when either 



e — >• or the number of nodes in the graph m — >• oo. In Section 2.4 we discuss F-convergence in more 
detail. Its most important feature is that if a sequence of functions {fn}^=i T-converges to a limit 
function /oo and in addition satisfies a specific compactness condition, then minimizers of /« converge 
to minimizers of /oo. 



It has been proven ga HHl |16|3 that ^ F-converges as e — )• to the total variation functional 

F^^{u) := a{W) [ \Vu\ (1.2) 



where now u is restricted to functions of bounded variation taking on two values (corresponding to the 
minima of the potential W) almost everywhere and the surface tension coefficient a{W) is determined 



by the potential W (see Section 5.1 for more details). Because the total variation of a binary function 
is proportional to the length of the boundary between the regions where the function takes on different 
values, from this limit functional the phase separating behavior can be seen clearly: u takes on one 
of two values, corresponding to the two different phases of the material and by minimizing the BV 
seminorm of u the interface between the two phases gets minimized. 

One of the results in this paper is a similar F-convergence statement for the graph Ginzburg-Landau 
functional: 

m ^ m 

feiu) ■=xYl ~ ^i)^ e ^ W{ui), (1.3) 

i,j=l i=l 

where Uj is the value of u on node i, ujij the weight of the edge connecting nodes i and j, m is the 
number of nodes in the graph, e > and x ^ (0) cc) is a constant independent of e and m, usually 



chosen to be x = | so the first summation is the analogue of / jVup (see Section 2.2). The different 
terms in this functional and its scaling will be explained below. 

The Euler-Lagrange equations for this functional are a nonlinear extension of the graph heat equa- 
tion using the graph Laplacian ^26j. Nonlinear elliptic equations on graphs were investigated in j48j 
and recently in j44J their well-posedness was studied. 

We study not only the limit e — t- in analogy with the classical continuum result, but also investigate 
the limit m — t- oo. For a graph embedded in M" this can be interpreted as the limit for finer discretization 
or sampling scale. In order to make sense of this limiting process we need to assume some additional 
structure on the graph that tells us how nodes are added along a sequence of increasing m. In this 
paper we consider 4-regular graphs {i.e., each node is connected to exactly 4 edges) with uniformly 
weighted edges in Sections |4] and [5j and a completely connected graph for the nonlocal means functional 
as studied in Section [6j but it is an interesting question if and how this can be extended to different 
types of graphs. Adaptation of our results to a 2-regular graph is fairly direct, but it is not clear at 
this moment how to extend our method to other graphs, even regular ones. 



^As poster child for F-convergence the proof has been reproduced, clarified, and extended upon in various ways, see 
e.g. [ISl[Sl[Mllll[nilll[31[lZllll 



r-CONVERGENCE OF GRAPH GiNZBURG-LANDAU FUNCTIONALS 



3 



1.2 Different scalings on a 4-regular graph 



The formulation of fg in (1.3) does not require the graph to be embedded in a surrounding space, 
although an embedding may exist as in the case of the 4-regular graph considered as an x square 
grid on the flat torus T^. 

We study two natural scalings for the functional on this 4-regular graph. The first is a direct 



reformulation of the graph functional fg from (1.3) with X = \ ^-iid weights equal to ^ on all 



existing edges and zero between two vertices that are not connected by an edge: 

N N 

hN,e{u) := N-^ ^ {ui+ij - Uijf + {uij+i - Uijf + e"^ ^ W{uij). (1.4) 

i,j=l i,j=l 

The second we get from discretizing the functional F^^ on the square grid using a forward finite 
difference scheme for the gradient and the trapezoidal rule for the integrals: 

N N 

kN,e{u) := e ^ - Uij f + (wjj+i - Uijf + e^'^N'^ ^ W{uij). (1.5) 

i,j=l i,j=l 

The subscripts in Uij denote the horizontal and vertical coordinates along the square grid. 

We will consider F-limits of these functionals when e — )• and N ^ oo sequentially. We also 
prove results in the case where we set e = iV~" for q > in a specified range and take — )• oo. 
Based on the F-convergence result in the continuum case we expect h^^g and k^^g to converge to total 
variation functionals. This intuition turns out to be correct, but with a twist: k^^s converges to the 
total variation functional fj,2 \ Vu\, but /i7v,e converges to the anisotropic total variation Jj2 \ux\ + \uy\- 
It picks up the directionality of the grid. Precise results are stated and proved in Sections [4] and [5j 
These results fit in very well with the research on T-convergence of discrete functionals to continuum 
functionals, as in e.g. [151 El E ESI El E El [23] . In fact, many of the techniques used in Section [5] are 
inspired by [5] specifically. 

We like to point out that there is also a substantial literature on the convergence of graph Laplacians 
and their eigenvalues and eigenvectors to continuum limits. See e.g. [iQl [361 [El SD [HI [571 113] and 
references therein. The techniques used and the kind of results obtained in those papers are quite 
different from ours, but in a certain sense our results can be seen as nonlinear extensions of the graph 
Laplacian case. 

In all cases we have to impose extra constraints on minimizers of the Ginzburg-Landau functional 
to avoid trivial minimizers. We show prove results showing that in most cases the addition of a mass 
constraint or the addition of an fidelity term to the functional is compatible with the T-convergence 
results. 



1.3 Asymptotic behavior of nonlocal means 

The functionals of nonlocal means type — or (anisotropic) nonlocal total variation type — we consider 
are built on the square grid in which the graphs are fully connected, [211 [Ml [351 [20] . Fix ^ G C°°(T^). 
We study 

N 

gN{u) := N~'^ ^ {uJL,N)i,j,k,i\ui,j - Uk,i\, (1.6) 

i,j,k,l=l 
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where lol,n ■= e~ l,n/'^ with a,L > constants (possibly depending on N) and 

r,s=—L 



If $ is thought of as an image on T^, as in e.g. [Uj, then L gives the size of the pixel neighborhoods 
whose pairwise comparisons form the graph weights. As we will see in Section|3j qn arises as the F-limit 



of /e in (1.3) as e — )• on this particular fully connected graph. This is a natural class of problems for 



which to study F-convergence as — t- oo. 



1.4 Structure of the paper 

This paper is structured as follows. Section [2] sets up notation and gives more background information 
about how to set up a PDE-to-graph 'dictionary' used to find the graph analogue of the Ginzburg- 
Landau functional. It also gives more details about F-convergence. In Section [3] the F-convergence 
result for is proved. This result holds for general finite undirected weighted graphs. Next we turn 
our attention to the square grid on the torus. In Sections |4] and [5] the F-convergence results for /ijv^ 
and kisf^s respectively are stated and proved. F-convergence for the nonlocal means type functional 
gjsf is discussed in Section [6j We close with a discussion of our results and open questions for future 
research in Section [71 



2 Setup 

We will start with introducing some general graph theoretical notation. 
2.1 Graph notation 

Let G = {V,E) be an undirected graph with vertex (or node) set V, \V\ = m £ N, and edge set 
E C V"^. Consider the space V of all functions 1/ — ?• M. A function u £ V can be seen as a labeling of 
the vertices of G. It is useful to number the vertices in V from 1 to m (in arbitrary but fixed order). 
We will write Im for the set of integers i satisfying 1 < i < m. IfuGV and rii £ V is the z*^ vertex 
we will use the shorthand notation Ui := u{ni). Let £ be the space of all functions — )• M, which 
are skew-symmetric with respect to edge direction, i.e., ii ip £ £ and eij := {ni,nj) € E is the edge 
between the i^^ and j*^ vertex in V we write (fij := ^{eij) and demand (pij = — 95jJj Since the graph 
is undirected we have eij € E <^ eji G E. When no confusion arises we will abuse notation slightly 
and consider eij = ej|^ In this paper we consider weighted graphs, which means we assume there is 
given a function uj : E ^ (0,oo), called the weight function, which assigns a positive weight to each 
edge. Because the graph is undirected the weight function is symmetric: ujij := uj{eij) = tOji. It is often 
useful to extend cj to a function on V'^ instead of on £^ C V'^ by identifying the edge eij with the pair 
(rij, Tij) G V'^ of its end vertices and setting uJij = if and only if eij ^ E. In particular, if the graph 
has no self-loops, ua = for all i £ 1^. In the same way we can extend (p £ £ to a function (/9 : — )• M 

^We impose skew-symmetry so that £ can be viewed as the space of flows as defined in e.g. 22] Section 2.2]. An 
interesting topological structure arises in this setting [221 Section 3] , but for our current purposes the demand of skew- 
symmetry neither hinders nor helps us. 

*See note|3] 
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by setting it to zero on node pairs that are not connected by and edge. We can incorporate unweighted 
graphs in this framework by viewing them as weighted graphs with the range of uj restricted to be 
{0,1}. We define the degree of vertex nj as di := X^jg/^Wjj. If G has no isolated vertices, then for 
every i £ Im di > 0. 

2.2 Graph Laplacians, Dirichlet energy, and total variation 

Our first goal is to define operators that serve as the graph gradient and graph divergence operators. 
Using these operators we can then define a graph Laplacian, a Dirichlet energy, and isotropic and 
anisotropic total variations on the graph. There are many possible choices to do this. Ours follow [4H 
Section 2] and [35] and are presented here. In Appendix |A] we give details and background on the 
justification of these choices. V = M"^ and[^£^ = M™('"~i)/2 are Hilbert spaces defined via the following 
inner products: 

{u,v)v ■■= ^ UiVidl, {(f, (/))£:= - J2 fij^ij^lr^^ 

for some r G [0, 1] and q € [1/2, 1]. Different choices of r and q are useful in different contexts, as will 
become clear later in this section. We also define the dot product as operator from £ x £ to V for 

If , (j) £ £ as 

{if ■ (p)i := 2 'Pij'l^ii^'ir^ ■ 

With the Hilbert space structure in place, if we define a difference operator, then all the other operators 
and functionals will follow naturally. We define the difference operator or gradient V : V — t- <f as 

(Vn)y := uj]p{uj - Ui). 

Notice that the choice q = 1 makes the gradient operator nonlocal on the graph, because its dependence 
on ujij disappears. The locality reappears in the if-'inner product' (or strictly speaking sesquilinear 
form), which is thus turned semi-definite. The opposite is the case for q= \. 
The other graph objects of interest for this paper now follow: 

• Norms: 

- ||u||v := ^/{u,u)v = ySi/T^S' 

- Ilv^lli :=y(^^^= VT^jeZ!^^!^^^^ • Note that II • ||£-,dot e V, 

- ||^i||v,oo : = max{|uj| :ie/m} and ||9j||£-,oo : = max{|99jj| -.ijGlm}- 

• The Dirichlet energy does not depend on r or q: 

^The factor | in m(m — l)/2 comes in because the graph is undirected. Strictly speaking £ is not isomorphic to 
l^m(m-i)/2 j£-^g ijjipose skew-symmetry, but this distinction is not relevant for our purposes. The Hilbert space structure 
can be defined in any case. 
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The divergence div : S ^ V defined as the adjoint of the gradient]^ 



• A family of graph Laplacians A,. := div o V : V — t- V (not to be confused with the p-Laplacians 
from the hterature). Writing out this definition gives 

(A,.M)i := d]-'ui - '^Uj = '^{ui - Uj). 

If we view u as a vector in we can also write 

ArU = {D^~'' - D-''W)u, 

where D'^ is the diagonal matrix with diagonal elements Da = d\ and W is the weight matrix with 
elements Wij = uJij. We can recover two of the most frequently used graph Laplacians from the 
literature (c/. |26l 1561 HT] ) by choosing either r = or r = 1. For r = we get the unnormalized graph 
Laplacian, for r = 1 we have the random walk Laplacian, which also goes by the name of (asymmetric) 
normalized Laplacian. For the latter case, the connection with random walks comes from the fact that 
D~^W is a stochastic matrix, i.e., the sum of the elements in each of its rows equals 1. Note that A,, 
is only symmetric if r = 0. A third graph Laplacian which is often encountered in the literature is the 
symmetric normalized Laplacian I — D^2WD^2^ where I is the m by m identity matrix. However, 
this one does not fit well into the current framework and we will not consider it here. 

• Total variation^ 

- The isotropic total variation TV : V — >■ M defined by 

TV(n) := max{(div(/?, ii)v : (p £ S, max WfWi < 1} 



A family of anisotropic total variations TVag : V — t- M defined by 

TVag(n) := max{ (div V : (f e£, \\f\\£,oo < 1} = ^ ^ -Uj\- 

TV and TVai appear in [35] as isotropic and anisotropic total variation respectively. In this 
paper we show that TVai is the F-limit of a sequence of Ginzburg-Landau type functionals 



(Theorem 3.1). 



If the graph has an isolated node i for which uiij = for all j and di = we interpret this definition as {div(fi)i = 0. 
'^An interesting question which falls outside the scope of this paper is in which respects, if any, the curvatures derived 
as 'derivatives' from these total variations resemble the continuum case curvature. 
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Because we can identify u G V with a vector u G M"* and all norms on M™ are equivalent, we can 
express convergence in any of these norms. For definiteness we choose a simple norm, not dependent on 
the degree function d: For a sequence {un}^=i C V and Uqo G V and corresponding vectors u„, Uoo G M™ 
we define 

Un — ^ "Uoo as n — )• oo iff \un — U00I2 — )• as n — )• 00, 
where | • |p, with p G N, is defined for u £ as 

\u\p := ( ^ 

subscript i labeling the elements of the vector. 

Where this does not lead to confusion, we will use the same notation u for both the function ti G V 
and the corresponding vector u G M"*. 



2.3 The functionals 

A standard choice of double well potential is W{s) = s^(s — 1)^. This is a representative example 
in the class of potentials for which our results hold. We always assume that W G C^(M), W > 0, 
and W{s) = iff s G {0, 1}. Different lemmas and theorems in this paper require different additional 
assumptions: 

(Wi) There exists two disjunct open intervals Iq and Ii containing and 1 respectively, constants cq, 
ci > 0, and a /3 > 0, such that 

< max{VF(s) : s G /} < m.m{W{s) : s G where I := IqU h and 

ys e io W{s) > co\s\'^ and Vs G /i > ci|s - l]''. (2.1) 



(W2) There exists a c > such that for large |,s| W{s) > c(s^ — 1). 

(W3) There exist ci, C2 > and p > such that for large \s\ ci\s\f < W{s) < C2\s\p. 

(W4) There exist 03,04 > and q> such that for large \s\ C3\s\'^ < W'{s) < C4\s\'^. 

{Wi) describes the behavior near the wells. It says that W is strictly bounded away from zero outside of 
neighborhoods of its wells and inside these neighborhoods W has a polynomial lower bound. We need 
it when we study the simultaneous scaling F-limit for hfj- and gives us explicit estimates of how quickly 
sequences of functions with bounded Ginzburg-Landau 'energy' approach the wells of the potential. 
Assumption {W2) is a coercivity condition that will help establish compactness in some situations 
(it could be replaced by any assumption that allows the conclusion that Jj2 W{u) is bounded from 
below by a function which is coercive in ||^t||L2(']r2)). (W3) with p > 2 is a condition needed to prove 
compactness in the classical Modica-Mortola F-convergence result for F^^ (see e.g. |54l Proposition 



3]). In addition we will use its lower bound to prove equi-coerciveness (Definition 2.2) of the functional 
k'^, which is defined below in (2.4). Finally we will need (VF4) to control the behavior of W in between 
grid points, when studying the simultaneous scaling F-limit for kf^. As is easily checked the standard 
example W{s) = s'^{s — 1)^ satisfies all the above assumptions for correctly chosen constants. 
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We frequently encounter binary functions in V and write 

V'' := {u G V : Vi G m G {0, 1}} . 



The graph Ginzburg-Landau functional : V 
energy: 



from (1.3) can be defined in terms of the Dirichlet 



^ m 

fe{u) = 2x\\Vu\\l + W{ui), with X e (0, oo). 

^ i=l 

Let be the two-dimensional flat unit torus. We construct a square grid with m = N"^ nodes 
Gn ■■= N-^Z^ nT^. Interpreting Gn graph we can use the notation from Section |2.l| with 

subscript A^, e.g. Vn are the vertices of Gn, Vn are the real- valued functions on Vn, {V^)n the 
binary ({0, l}-valued) functions on Vj\[, etc. We understand the vertices Vn to be embedded in T^. To 
distinguish the horizontal and vertical directions in our graph when working on Gn, instead of Ui we 
will write Uij := u{nij) where riij := {i/N,j/N) G Vn C T^. We refer to a single square in the grid 
by ' ' 

:= [i/N, {i + 1)/N) X [j/N, (i + 1)/N). (2.2) 

We remind the reader that we introduced three different functionals on the square grid: the graph 
theoretical Ginzburg-Landau functional /itv.e : Vat — )• M in (1.4), the discretized Ginzburg-Landau 
functional /ctv.e : Vat — )• M in (1.5), and the 'sharp interface' {i.e., e — )• 0) nonlocal means functional 
gN ■■V^^^Rm^. 

We call hN,e the graph theoretical Ginzburg-Landau functional because it is equal to fs from (1.3) 
if we choose the weight u in as 



N'^ if {\i-k\ = l A j = l)\/{i = k A \j -l\ = 1), 
otherwise, 



and X = 5- f^N,e we get by using the trapezoidal rule and a standard finite difference scheme to 
discretize the Ginzburg-Landau functional F^^. We have used the periodicity to relate the terms of 
the form {uij —Ui-ij)'^ and {uij — Uij-i)'^ to {ui+ij — Uij)'^ and (ujj+i — Ujj)^ in the sum, respectively. 

To study L-convergence for the 'simultaneous' limits e — )• and — )• oo of hN,e and kN,e we set 
e = N~°', for a > 0, and let — t- oo in the functionals 



TV 



N 



h%{u) := A ^ ^ {ui+ij - Uijf + {uij+i 
N 



k%{u) 



A" 



"I- J 



N 



(2.3) 
(2.4) 



We prove L-convergence results for h'^ in Section 4.2 and for k'^ in Section 5.2 

Note that h'^ = N'^~^k]^ if 7 = This shows that we do not expect /i^ and k]^ to have the 

same limit, unless possibly if a = 7 = 1. This value falls outside the regimes for a we consider and 
hence we do find different limits. 
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2.4 F-convergence 

T-convergence was introduced by De Giorgi and Franzoni in [29j. It is a type of convergence for 
function(al)s that is tailored to the needs of minimization problems as we will see below. A good 
introduction to the subject is [l5], the standard reference work is [28] , 



Definition 2.1. Let X he a metric space and let {-F„}^-^ he a sequence of junctionals Fn : X ^ 

M U {±00} . We say that Fn T -converges to the functional F : X — t- M U {±00}, denoted hy Fj — t- F if, 
for all u £ X we have that 

(LB) for every sequence such that Un ^ u it holds that F(u) < liminf F„(n„,) and 

(UB) there exists a sequence {n„}^;^ such that F{u) > lim sup i^„(n„). 

The lower bound condition (LB) tells us that the values along the sequence F„(n„) are bounded 
from below by F{u), the upper bound (UB) shows that the value F(u) is actually achieved. Combined 



with a compactness or equi-coerciveness condition (Definition 2.2 below), this allows for conclusions 
on the minimizers of Fn and F. 

It is useful to note that to prove (LB) for a given sequence {un}^=i we only need to prove it for 
a subsequence {iin'}^=i C {u„}^^ such that lim Fn/{un') = liminf F„(u„). If (LB) is satisfied for 

such a sequence (which always exists), then F{u) < lim inf F„/ (ti„' ) = lim Fn'^u'^) = liminf i^(u„). 

n^oo 

Hence, when proving (LB) we will assume without loss of generality that {unj^i is such a sequence. 
The uniqueness of the limit then implies that it suffices to prove (LB) for any subsequence. Clearly we 
can also assume that liminf F„(n„) < 00 and hence it suffices to prove (LB) for a specific subsequence 
{'"n"}^/=i for which there is a C > such that («„//) < C. In practice this means that to prove 
(LB) we can assume a uniform bound on Fn{un) and freely pass to subsequences when needed. 

If we are working with functionals that depend on a continuous parameter, e.g. e — )■ or — t- 00, 
we have to prove (LB) and (UB) for an arbitrary sequence {en}'^=i with e„ — as n — 00 (or {Nn} 
with Nr, — )• 00 as n — )• 00). 



00 

n=l 



Definition 2.2. Let X he a metric space and let {Fn}^^-^ he a sequence of functionals Fn : X ^ 
MU {±00}. We say the sequence is equi-coercive if for every t S M there exists a compact set Kf C X 
such that for every n G N {u £ X : Fn{u) < t} C Kt. 

In practice equi-coerciveness is proved by showing that any sequence {n„}^^ for which Fn{un) is 
uniformly bounded has a convergent subsequence. 

F-convergence combined with equi-coerciveness allows us to conclude the following result. 

Theorem 2.3 (Chapter 7 in [28j and Theorem 1.21 in p^). Let X he a metric space, he a 

sequence of equi-coercive functionals : A — )• M U {±00}, and let F he the T -limit of Fn for n — )• 00. 
Then there exists a minimizer of F in X and mm{F{u) : u € X} = lim mf{Fn{u) : u G X}. 

n—^oo 

Furthermore, if {n„}^-^ C X is a precompact sequence such that 

lim Fn{un) = lim inf{F„(n) : u G X}, 

n— >oo n— >oo 

then every cluster point of this sequence is a minimizer of F . 
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Following [5] for our purposes it turns out it is often more useful to reformulate F-convergence in 
terms of the F-lower limit 

F'{u) := infjliminf i^„(un) : n„ — )■ n} 

n— >oo 

and the F-upper limit 

F"{u) := infjlimsup F„(n„) : n„ — t- n}, 

n—^OD 

|28l Definition 4.1]. It can be shown, [28, Remark 4.2, Proposition 8.1], [LQl, that our definition of 
F-convergence above is equivalent to the following two conditions. For each u £ X we have that 

(LB') F{u) < F'{u) and 

(UB') F{u) > F"{u). 

The benefit of this reformulation is that the functions F' and F" are lower semicontinuous \28\ Propo- 
sition 6.8], which comes in handy in Sectionjsj In fact, since conditions (LB) and (LB') are equivalent, 
we will sometimes use the combination (LB)+(UB') to prove F-convergence in this paper. Note that 
(UB) implies (UB'). 

2.5 Constraints 

It is common in semi-supervised learning applications to have a mass constraint or an additional term 
in the functional corresponding to a fit to the known data. Moreover such constraints are typically 
necessary to obtain nontrivial minimizers. We need to check that these constraints are compatible with 
the convergence. 

First consider the case of an adding a fidelity term of the form X\u — f\p to the functional, where 
/ G V is a given function (usually representing some known data to which the minimizer should be 
similar) defined on some or all of the vertices in V and A > is a parameter. If p agrees with the 
topology of the F-convergence we can use the property that F-limits are stable under addition of a 
continuous term or a sequence of continuously convergent terms [28, Definition 4.7, Propositions 6.20- 
21] to conclude that the Ginzburg-Landau functionals plus fidelity term again F-converge. We will 
summarize the results that are relevant for us in the following definition and lemma, based on the cited 
definition and propositions in [28j . 

Definition 2.4. Let X be a metric space and let {Fn}'^=i he a sequence of functionals Fn : X ^ 
M U {±00} . We say the sequence is continuously convergent to a function F : X ^ ML) {±00} if for 
every n G X and for every ij > there is an N G N and a 6 > such that for all n > N and v G X 
with \\u — v\\ < 6 we have \Fn{v) — F{u)\ < rj. 

Lemma 2.5. Let X be a metric space and let {Fn}^:^i sequence of functionals Fn ■ X ^ MUjitcxD} 
which T -converges to F : X ^ {±00} and {Hn}^=i a sequence of functionals Hn ■ X — MU {±00} 
which is continously convergent to : X — >• M U {±cx)}, then F^ + T-converges to F + H . If 
G : X — 7- M U {±00} is continuous functional, then F^ + G T-converges to F + G. 

If instead a mass constraint is imposed on the minimizers we have to check that each convergent 
sequence preserves the constraint (to make it compatible with the lower bound and compactness con- 
ditions) and that the recovery sequence from the upper bound condition does satisfy the constraint 
(or can be adapted to satisfy it, without violating the upper bound condition). Since we are dealing 
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with LP convergence, the former is usuahy triviahy satisfied, but the latter does demand some more 
attention. Details for each of the functionals are provided in the relevant sections. 

For functions n G V'' we carefully need to determine the form of our mass constraint. The constraint 
X^I^i Ui = M leads to shrinking support when m — )• oo which is unwanted. An alternative condition is 
an average mass constraint of the form ^ X^I^i ~ Note that X^I^i "^i take on only integer 
values between and ni and hence mM should be of that form as well, in order for the average mass 
constraint not to lead to an empty set of admissible minimizers. If we impose this for all m this is 
only possible if M = or M = 1, however specific subsequences of m can be able to satisfy this 
condition for different values of M {e.g. if Af = ^ and we consider even m). Hence the choice of M 
can constrain the subsequences of m which are admissible. In order to avoid the possible difficulties 
with the average mass equality, one can also impose an average mass inequality. Since the arguments 
for the mass equality easily generalize to an inequality, we will not discuss this situation further. 

3 F-convergence of 

3.1 F-convergence and compactness 



In this section we prove F-convergence and compactness for the functional : V — )■ M from (1.3). 



Theorem 3.1 (F-convergence). /e — >• /o as e ^ 0, where 

f / xEije/^ ^ijWi - «il e V^ ^ r 2xTVai{u) tf u e V^ 

° I -|-oo otherwise \ +oo otherwise. 

Theorem 3.2 (Compactness). Let W satisfy the coercivity condition {W2), let C M+ be a 

sequence such that En ^ as n ^ 00, and let {un\'^=i C V be a sequence such that there exists a 
C > such that for all n £ N fe„{un) < C. Then there exists a subsequence {wn'}^=i C {unjj^i and 
a G such that Un' — >• ^^oo as n —t- 00. 

Although in F^^ the first term is scaled by e, the first term of /e contains no e. The reason for this 
is that the Dirichlet energy in Ff^ is unbounded for the binary functions u that form the domain of 
the limit functional Fq^. However, the difference terms in are finite even for the binary functions 
and thus need no rescaling. The proof of F-convergence uses this fact, to view the difference terms as 
a continuous perturbation of the functional 

We{u) := ^ X] ^(^*)- 

Lemma 3.3. The sequence of functionals T-converges: — )• wq where wq :V ^ {0, 00} is defined 

e-s>0 

via 

[ -|-oo otherwise. 

Proof. To prove the required lower bound (LB) let u G V and consider sequences and {un}'^=i 

such that £n —5- and — )■ n as n — )• 00. If u G from g^ > {) \i follows that wq{u) = < 
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liminf We„{un)- If on the other hand u G V\V then for large enough n 3v such that Un{v) ^ {0, 1} 

n— >oo 

and hence ^ 

liminf W£^{un) > hminf — W{v) = oo = ■wo{u). 

n— ^oo n— >oo Sn 

For the upper bound (UB) we can assume without loss of generality that n G V^. Define for every 
n G N Un ■= u, then trivially n„ — )• « if n ^ oo and furthermore limsuptt;e„(n„) = = wo{u). □ 

n— >oo 

Proof of Theorem 13. li Define 

w{u):=x ^ uJij{ui-Ujf. 
w is a polynomial on M™ and hence continuous. F-convergence is stable under continuous perturbations 



{e.g. |28l Proposition 6.21]). Since fs is a continuous perturbation of We we have by Lemma 3.3 
/g -!-> u) + ifo as e — )• 0. We complete the proof by noting that if n G V'', then 

W{u) = X ^ UJijlUi -Uj\. □ 

Proof of Theorem 13.21 By the uniform bound on ff,^{un) we have 



W[{Un)i) < Cen. 



Combined with the coercivity condition {W2) on W we conclude that for n large enough the sequence 
{wnj^i is bounded and hence by the Bolzano- Weierstrass theorem there exists a converging subse- 
quence with limit Uoo- Since W G C^(M) and e„ — )■ as n — )• 00 we conclude that W {uoo{vi)) = for 
all i G Im and hence Uqo G V^. □ 

Remark 3.4. Since /o is defined on binary functions n, if we write 

Sk ■■= {i £ Im ■ Ui = k} for k G {0, 1}, 

we can rewrite /o as 

[ -|-oo otherwise 

Minimizing /o thus corresponds to finding a minimal graph cut of G, i.e., dividing the graph into 
clusters with minimal edge weight between them. Such a minimization requires an extra constraint to 
avoid trivial minimizers. A common choice is to prescribe the number of clusters (in this case two) 
one wants, or to introduce a normalization into the sum of the weights based on the cluster sizes {e.g. 
normalized cut, normalized association [53], and Cheeger cut |55j)). One could also minimize /o under 
a fixed mass constraint (see Section 3.2 below). One often relaxes the problem of normalized graph cut 
minimization by losing the binarity constraint, e.g. in spectral clustering |49|. 156]. 

If the weight function ui is such that there is a nontrivial partition Au B = Im such that Uij = 
for i £ A and j £ B then a nontrivial minimizer of /o (without the mass constraint) is clearly given 
via 5*0 = ^ and Si = B. In this case we can write 

Ui — 

di 

which is of a form related to the image denoising method known as nonlocal means, of. 



r-CONVERGENCE OF GRAPH GiNZBURG-LANDAU FUNCTIONALS 



13 



Remark 3.5. In this paper we assume that W has wehs at and 1. The proofs in this section 
make no use of this fact and can easily be extended to potentials W with wells at values si and S2- In 
this case the set needs to be redefined as the set of functions taking values in {si, S2} and the limit 



functional /o from Theorem 
for u G V^. 



3.1 



is multiplied by a factor |si — S2I, since {ui 



\Sl - S2\\Ui 



Ui 



3.2 Constraints 



Next we show that the addition of a fidelity term X\u — f\p or a mass constraint is compatible with the 
F-convergence. 

Theorem 3.6 (Constraints). + A| • —f\p — )• /o + A| • —f\p as e — >• 0, where p G M+, A > 0, and a given 
function f £ V (or possibly a given function / : [/ — t- M where U is a strict subset of the vertex set V 
and the sum in \ ■ —f\p is restricted to vertices in U ). Compactness for + A| • —f\p as in Theorem 3.2 
holds. If instead, for fixed M > 0, the domain of definition of fe is restricted to 



V 



M 



{u G V : ^ Ui 



mM}, 



where M is such that mM is an integer between and m, then the results of Theorems 3.1 and 3.2 
remain valid, with the domain of fo restricted to V^^ . 

Proof. The fidelity term A|u — f\p is a polynomial, hence a continuous perturbation independent of e 



to fe and thus F-convergence follows by Theorem 3.1 and Lemma 2.5, The addition of this term does 
not affect the compactness property at all. 

The mass constraint is compatible with the limit functional being defined on binary functions. The 
constraint is preserved under convergence in V = M™, so it is compatible with (LB) from Definition 2.1 



and compactness. If G V satisfies the mass constraint, then trivially so does the recovery sequence 
{wn}^! used to prove (UB) from Definition 2.1 in Lemma 3.3 □ 



4 F-limits for the graph based functional h 



In this section we will study the convergence properties of /ijv,e from (1.4). We consider two different 
cases. In the first we first take the limit e — )• and then N — )■ 00, in the second we take both limits at 
once by substituting e = N^"' for well chosen a > and then considering the limit N ^ 00 for hfr in 



(2.3). These results have a similar feel as numerical convergence results, however the lack of regularity 



of the binary limit functions complicates the results and proofs. 



Remark 4.1. We note that, while we give the proofs for T in Sections 4.1 and 4.2 they can easily 



be generalized to T for any d G N, if we let the scaling factor in the first term of hjsf^g in (1.4) be N 



instead oi N ^ and we change h^ in (2.3) accordingly 



rl-d 



N N 

hN,e{u) := A^^""^ {ui+i,j - Ui,jf + {uij+i - Uijf + e"^ ^ W{ui,j), 

i,j=l i,j=i 
N N 

h%{u) := Afi-^ Yl - ^^.i)' + - + N'' ^ W{ui ' 

i,j=l i,j=i 
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For the extra fidehty term in Theorem 4.13 the scaling factor then needs to be instead of N 



4.1 Compactness and sequential F-limits: first e 0, then N oo 



p 

By Theorem 3.1 we immediately have /iat^e — ?• /iat^o as e — )• 0, where /iat^q is defined for u G V as 



hN,o{u) :-- 



N ^ Eij=i - Uij\ + \uij+i - Uij\) if M G V^, 

+00 otherwise. 



Jj2\u^\ + \uy\ if^xG W(T2;{0,1}), 
+00 otherwise, 



In this section we prove that /iat^q — ^ /ioo,o as — )• oo where /ioo,o is defined for u G L^(T^) as 
hoo,o{u) := 

The anisotropic total variation in /ioo,o is defined as 

/ |na;| + |uj^| := sup < / udivv : V £ CI{T'^;M.'^),\/x \v{x)\nr, < 1 

where for a vector v{x) = {vi{x) , V2{x)) G M^, the norm | • |oo is defined by 

\v{x)\oc ■■= max{|fi(x)|, |v2(x)|}. 

For functions u G BV{T'^; {0,1}) this anisotropic total variation gives the length of the (reducecQ 
boundary of the set {u = 1} projected onto the horizontal and vertical axes (counting multiplicities). 
We prove a compactness and F-convergence result. 

Theorem 4.2 (Compactness). Let {Nn}'^=i C N satisfy Nn ^ oo as n ^ oo and let {un}'^=i C L^(T^) 
be a sequence for which there is a constant C > such that for all n G N hN„,o{un) < C. Then there 
exists a subsequence {'Un'}^=i C {n„}^]^ and a u £ BV{T'^; {0,1}) such that Un' u in L^(T^) as 
n' — )• oo. 

Theorem 4.3 (F-convergence). hj^fl — )■ /ioo,o as N ^ oo in the L^(T2) topology. 

The convergence of /iat^q to an anisotropic total variation is reminiscent of the related, but different, 
results in [23] . 

Because the limit function /ioo,o is defined on L^(T^) functions, it will be useful to identify the 
binary functions on the square graph, i.e., the functions in V^, with a subset of L^(T^), namely binary 
functions on that are piecewise constant on the squares of the grid. Using the notation S]^ from 



(2.2) we define 

An '■= G L^(T^) : y{i,j) £ 1% u is constant a.e. on | , (4-1) 
An := {uGAn -.uG L^(T^; {0, 1})} . 

We construct a bijection between the normed spaces Vat and by identifying u G Vat with the unique 
u G An which satisfies u = Uij a.e. on S^j^ for all {i,j) G It is easy to check that convergence 

*For the definition of reduced boundary see [3 Definition 3.54]. 
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in Vn corresponds to convergence in An (1 < P < oo) and that the bijection maps the subset 
to A% and vice versa (for its inverse) . In what foUows we will drop the tilde if this does not lead to 
confusion. 

With this identification we write for n G 

(u) I '^'^^ + 1^1/1 if "U G -^Af' 2) 

I +00 otherwise. 

In fact for our F-convergence purposes without loss of generality we extend the functional to all u G 
L1(T2), such that hNfi{u) = +oo if u G ^^(T^) \ A%. 
First we prove the compactness result. 



Proof of Theorem 14.21 By the definition of the isotropic and anisotropic total variation and (4.2) 
we have for all n G N 



/ |Vti„| < / \{Un)x\ + \{Un)y\ < C. 
JT2 JT2 



In addition for each n G N G A^j^^, hence 11^^11^1(12) < 1- We deduce that the sequence {u„}^^ is 
uniformly bounded in the BV norm and thus by compactness ( j37l Theorem 1.19] or [62\ 5.2.3 Theorem 
4]) there exists a subsequence {Mn'}Jv=i {un}^=i and a li G BV{T'^) such that Un' — )• u in L^(T^) as 
n' — )• oo. Since all Un take the values and 1 a.e. by pointwise a.e. convergence (after possibly going 
to another subsequence) so does u. □ 



In the next lemma (LB) from Definition 2.1 is proved. 



Lemma 4.4 (Lower bound). Let u G ^^(T^) and let {un}n=i C ^^(T^) and {iVn}^=i C N 6e such that 
Un ^ u in L^(T'^) and Nn — )■ cxd as n ^ oo. Then 

/ioo,o(^i) < limmfhN^ oiun). 

Proof. First consider the case where u G BV{T'^; {0, 1}), then without loss of generality we can assume 
that Un G A'j^^ . Analogous to the isotropic total variation also this anisotropic total variation is lower 
semicontinuous with respect to convergence |24l Lemma A. 5] hence we find 

hoofiiu) = / \ux\ + \uy\ < lim inf / |(n„)a;| + |(u„)j^| = liminf /iAr„ o(^tn)- 

If n G L\T^) \ BV{T^; {0, 1}) and Un ^ u in L\T^), then 

hoofl{u) = oo = lim inf hN^fiiun), 
which we prove via contradiction: Assume lim inf /i7v„,o(^n) < oo, then there is a subsequence C 



{un}'^=i for which hN^,fi{uni) is uniformly bounded and hence by Theorem 4.2 — )• u in L^(T^) as 



n' — )• oo, where u G BViT ] {0, 1}). By the uniqueness of limit u = u which is a contradiction. □ 

To prove the lim sup inequality we use the following results from [2l] (which we give here in a form 
adapted to our situation which follows easily from the results in 
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Lemma 4.5 (Corollary A. 4 and Theorem 4.1 from [24J). For u G BViT"^) we have 

I '^X I ~t~ I Uy I 

sup I j udivv : v £ L'^{T'^;R'^), divv £ L'^{T'^), \v{x)\oo < la.e.}. 
Furthermore, for each u G A^j^ there exists a v £ L°°(T^;]R^) such that \v{x)\oo < 1 a.e., 



u divv = / \ux\ + \uy\, 

T2 



Il°°(t2) — 4iV. 



The first result in the above lemma says we can relax the condition on the admissible vector fields in 
the definition of the anisotropic total variation to L°° vector fields with essentially bounded divergence. 



The second result shows that the supremum is achieved by a specific vector field if u G Ji'j. 



Lemma 4.6 (Upper bound). Let u G -L^(T^) and {Nn]'^^^ C N 6e such that Nn — t- cxd as n ^ oo. 
Then there exists a sequence {un}'^^i C -L^(T^) such that Un ^ u in L^(T^) as n —t- oo and hoofl{u) > 
lim sup hNr^oiun)- 

n—^oo 

Proof. Without loss of generality we can assume that u G BV(T'^;{0,1}). Construct Un as follows. 
For X G 5]y^ define 

^^(^) — I 1 if (SnJ^ C supp-u, 
" \ otherwise. 

In words, Un takes the value one on those squares of the grid whose interior lies completely in the 
support of u and zero on the other squares. By 9=i,suppu denote the reduced boundary of the set 
suppn, i.e., all points in Ssuppu for which there is a well defined normal vector (see [9, Definition 
3.54]). Since suppu^ is the maximal set which is both a union of squares on the grid and is contained 
in supp u, the difference in area between supp u and supp n„ is bounded by the length of the reduced 
boundary of suppn times the area of a square, i.e., 



/ \un - u\<n^{d* supp u)N~^. 
Jt2 



Because u G ^^(T^; {0, 1}) the set suppu has finite perimeter and hence n„ — )• n in L^(T^). 

For each n„, let f„ G L°°(T'^;M?) be the vector field whose existence is guaranteed by Lemma 4.5 
then 



\iun)x\ + \{un)y\ = - n„ div -[;„ = - / udivvn+ {u-un)d[vvn- (4.3) 
t2 Jt^ Jt2 Jt2 



For the last term we have 



(n - Un) divvn 



< ||n - u„||ii(ir2)||divfri||L°o(T2) 

< An^id* suppu)iV^^ ^ 0, as n — 7- oo. 
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Hence, by the first statement of Lemma 4.5 we liave 



limsup / + |(^^n)?/| = limsupf — / ndiv 

< limsup 



Ux\ + = / \Ux\ + \Uy\, 



which proves the result. 



□ 



Proof of Theorem 4.3 Combining Lemmas 4.4 and 4.6 proves the F-convergence result in Theo- 
rem Ol □ 



Remark 4.7. Note that we could have used any space instead of in the results above. Because 
is bounded, convergence in implies convergence in and so the result of Lemma 4.4 is easily 



recovered. For Theorem 4.2 and Lemma 4.6 we note that because u„ and u are binary functions taking 
values and 1 a.e., the bound on their difference is the same as that on their difference and the 
results follow again. 

We end this section with an illustration of the preference for squares and rectangles of hjsffi. 

Lemma 4.8 (Minimizers of hNfl)- Let M G [0, 1] be such that N^M = for some K £ N. 

// M E [0, then uq is a minimizer of h^fl over all u G .4^ that satisfy fj2 u = M if and only if 
uq is the characteristic function of a square. 

IfM e (|,1) then Uq is a minimizer if and only if it is the characteristic function of a rectangle of 
the form R = [a, h\ x [0, 1] C or R = [0, 1] x [a, b] C for a, b that satisfy the mass constraint. 

If M = ^, Uq is a minimizer if and only if it is the characteristic function of a square or a rectangle 
R as above. 

Proof. First consider the square grid Gat to be embedded in instead of in so that we do not 
have periodic boundary conditions on [0, 1]^. Let u be the characteristic function of a set $7. We can 
assume is connected, because else hiyfl{u) can be lowered by rearranging Q, to be connected without 
changing the mass. Let uq have the square with sides of length Lq = KN~^ = \fM as support and let 
O be contained in a rectangle with sides of lengths L and B. Then YYn< LB and hence 

/i7V,o(^io) = 4Lo = 4(^y"^ n) ^''^ < aVlB < 2{L + B) < hN,o{u), 

with equality if and only \i L = B = Lq. Hence characteristic functions of squares are the minimizers 
of hNfi if we ignore periodic boundary conditions. 

However, on the periodic torus if is a rectangle we can use periodicity to eliminate two sides of 
the rectangle. This can be done only if the other two sides have length 1 and hence hN,o{u) = 2. This 
beats the square if 4Lo = 4\/M > 2. □ 



It is worth noting here that this asymptotic behavior of h^^^ is not accessible via numerical sim- 
ulations of a gradient flow, since it is dependent on e being small enough for the minimizers to be 
essentially binary and hence not differentiable. 
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4.2 Simultaneous scaling F-limit for h'^ 



N 



In Section 4.1 we first took the limit e — t- for /i^v.e before letting N ^ cc. In this section we will 



consider the limit if we let both parameters go to their limit simultaneously. To this end we choose 
£ = N~°' for a > and consider the limit N 
sense of Section 



oo of hfr in (2.3). We identify Vat with An in the 



5.1 



and extend hf^ to all of ^^(T^) by setting h%{u) := +oo for u G L'^{T^) \ An- 
We show that the F-limit is again given by /ioo,o if is large enough, depending on the growth rate 
(3 of W around its wells given in assumption (VFi). 



Theorem 4.9 (F-convergence) . Assume that W satisfies condition {Wi) for some (3 > in (2.1) and 
let a > (3, then hf^ — )• /ioo,o as N ^ oo in the L^(T^) topology. 



Theorem 4.10 (Compactness). Assume that W satisfies condition (Wi) for some (3 > in (2.1) and 

let a > (3. Let {Nn}^=i C N satisfy Nn — )• oo as n —)• oo and let {un}'^^i C L^(T^) be a sequence for 
which there is a constant C > such that for all n £N hfj < C. Then there exists a subsequence 
{un'}^=i C {un}^=i and au£ BV{T^; {0, 1}) such that Un' u in L^(T^) as n' — t- oo. 

We first prove the lower bound for the F-limit. 

Lemma 4.11 (Lower bound). Let u G L'^{T'^) and let {un}'^=i C ^^(T^) and {iVn}^=i C N be such 
that Un ^ u in L^(T^) and Nn — t- oo as n ^ oo. Assume that W satisfies condition (Wi) for some 
(3 > in (2.1), and let a > (3. Then 

hoofiiu) < liminf/i^ (n„). 

n^oa 

Proof. First consider the case where u E BV{T'^; {0, 1}). Without loss of generality we may assume 
that h'^^ is uniformly bounded. Since W is nonnegative we deduce there is a C > such that for all 
{i,j) G W{{un)ij) < CN~°^. Together with (2.1) in assumption {W2) this implies that for n large 

enough and all (i,j) G we have {un)ij G /. In addition, by the growth condition in (2.1) we get 
for n large enough and s G {0, 1} 



Cs\{u. 



f < Wi{u. 



"njt,j , 



< CN- 



Hence, if we define Sn '■= iC j min({co, ci})) A'^n we deduce that for n large enough bn G (0, 5) and 
for aU (i, j) G (urt)j,j G (-(^„, (5rt) U (1 - 1 + (Jrt)- Define 

■= {(hj) G In„ ■ ((^in)j,i G (-(5„,(5„) A (n„)i+ij G (1 - 5„, 1 + 5,,)) 
V ((u„)i+ij G ( ) A(u„)ij G (l-(5n,l + (5n))}, 

then for (z,j) G X„ we have \{un)i+ij - iun)i,j\ 



(4.4) 



> 1 — 25n, hence 



Y] (K)i+1„ 



{'^n)i,j) 



(4.5) 



>N~\l - 26n) \{Un)i+l,j - Mj\ + Yl 

N„ 



{{U'n)i+l,j - {Un)i,jY 
n)i+l,j ('Un)j,j| 
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For the second summation in ( |4.5[ ) we note that there are c > 0, C > 0, such that 

< C{N^-"/l^ + iV^"^°/^) ^0 as n ^ oo. 



The third and fourth summation in (4.5) we combine into 

^n^ \iUn)i+l,j - {Un)ij\{\iUn)i+lJ - iUn)i,j\ - 1) 

<N~^ ^ \{Un)i+l,j - {Un)i,j\ \\{Un)i+l,j - {Un)ij\ - l\ 

< N-^Nl26n{25n - 1) < CiN}-^""'^ - iV^~"/^) ^0 as n 



oo, 



for some C > 0. As in Section 5.1 we now identify Vn^ with Am^ to write for the first summation in 



^n^ Y - {Un)i,j\ 



/ l(^n):J 



hence from (4.5) and the computations that fohowed we deduce 



(4.6) 



Analogously we have 



N, 



N- 



-^"71 p 



By the lower semicontinuity of the anisotropic total variation with respect to convergence 
Lemma A. 5] we have 

/ In^'l + < liminf / \{un)x\ + \{ur 

Jj.2 n-s>oo Jj2 



^njy\i 



hence 



/ \ux\ + \uy\ <liminfiV„ ^ ((n„)i+i j 
< liminf h%^{un), 



{Un)i,jf + {{Un)i,j+1 - {Un)i,jy 



which proves the lower bound for the case where u G BV{T'^; {0, 1}). 

Now consider u G L^(T^) \ BV{T^; {0, 1}). Let m„ — )■ u in L^{T'^), then an argument from contra- 
diction (similar to that at the end of the proof of Lemma |4.4[ ) and the compactness proved below in 
Theorem |4.10 prove that /iqo o{u) = oo = liminf («„)• □ 
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Next we prove the upper bound. 



Lemma 4.12 (Upper bound). Let u G L^(T^) and {Nn]'^=i C N 6e such that Nn — >• oo as n ^ oo. 
Then there exists a sequence {un}^=i C L^(T^) such that Un ^ u in L^(T^) as n ^ oo and 

hoo,o{u) > limsup/i^^(u„). 

71— >00 

Proof. Without loss of generahty we can assume th at u G i?F(T^; {0, 1}). Construct Un G -^Nn as 
follows. For X £ S^j^^ define, as in the proof of Lemma 



4.6 



Unix) :-- 



1 if C supp-u, 

otherwise. 



Then {{un)i+i,j - {un)i,j) = \ {un)i+i,j - {un)i,j\ and hence identifying Vat^ with An„ 

\{Un)x\- 



Nn Nn 

2Z - {'^n)i,j) = ^ \{Un)i+l,j - ( 

i,j=l i,j=l 



U: 



'njt,j I 



f2 



Similarly 



Nn 



^ ^ {{Un)i,j+l 
«J=1 



{Un) 



\{Un)y\ 



Since every Un takes values in {0, 1} we can repeat the argument from the proof of Lemma 4.6 in and 



following (4.3) to prove that 



Since for every n G N and G , W{{un)ij 



limsup / \{Un)x\ + \{Un)y\ < / | | + | Wj/ 1 • 

n-s>oo Jj2 Jf2 

we get the desired result. 



□ 



Proof of Theorem 4.9| , Combining Lemmas 4.11 and 4.12 we get the P-convergence result in Theo- 
rem |49l □ 
Next we prove compactness. 

Proof of Theorem 4.10[ By the first part of the proof of Lemma 4.11 we have, after possibly going 
to a subsequence, that for all n G N and all (i, j) G {'^n)i,j S I, hence ||ti„||j;^i(']p2) is uniformly 

bounded. By the same proof, in particular (4.6) and the uniform bound on we have for all 

n G N that fj.2 \{un)x\ + |(^n)j/| is uniformly bounded. We deduce as in the proof of Theorem |4.2| a 
uniform bound on the BV norms of Un from which it follows by the compactness theorem in BV f|37t 
Theorem 1.19] or [321 5.2.3 Theorem 4]) that there exists a subsequence {^in'}^=i C {un}^=i a nd a 
u G BVir"^) such that Un' — )• u in L^(T^) as n' — )• oo. By the arguments in the proof of Lemma 4.11 



each Un' takes values in (— (5„',(5„/) U (1 — 6n',^ + 5^') where (5„/ — )■ as n' — oo. After possibly 
going to another subsequence Un' converges pointwise a.e. to u, hence u takes values in {0, 1} almost 
everywhere. □ 



r-CONVERGENCE OF GRAPH GiNZBURG-LANDAU FUNCTIONALS 



21 



4.3 Constraints 

In this section we show that addition of a fidehty term to the functional or imposing a mass constraint 
are compatible with the three F- limits we discussed, i.e., e — t- for h^^e, N ^ oo for h^fi, and N ^ oo 
for h%. 



Theorem 4.13 (Constraints). (1) hN,e + AA^"^| • -f\^ 4 hN,o + AA^^^j . _^|p ^ ^ g, where pGN, 
A > 0, and a given function f € Vat (or possibly a given function / : [/ — t- M where U is a strict subset 
of the vertex set V and the sum in \ ■ —f\p is restricted to vertices in U). A compactness result for 
+ XN-^\ 



h 



-f\p as in Theorem 



3.2 



holds. 



If instead, for fixed M > 0, the domain of definition of h^^e is restricted to {i.e., V^^ from 



Theorem 3.6 on the grid Gjy) where M is such that N'^M is an integer between and N"^ , then the 
T -convergence and compactness results for e — )• remain valid, with the domain of h]\jfl restricted to 
V*^ as well. 

(2) Let p £ N, X > 0, f £ C^(T^) and f^ G An the sampli ng of f on the grid Gjv (/, /iv 0''nd 
their norms can also be defined on subsets o/T^ and Gn as in part 

^00,0 + A|| • —f\\\p^j2-^ as N ^ oo in the LP(T^) topology. A compactness result as in Theorems 
holds for hj\fQ + AiV^^| 



Remark 



4.7 



TV 



4.2 



and 



-Jn\p- 



If instead the domain of h]\!fi is restricted to , for a fixed M such that N M is an integer 

and 



between and N"^ , then the compactness and lower bound results from Lemma 4.4, Theorem 4.2 
Remark 4.7 remain valid, with the domain of /ioo,o restricted to 



Wm(T2; {0, 1}) := [u G ^(T^; {0, 1}) : J^^u = m}. 



With the same restriction, the upper bound result from Lemma 4.6 is still valid if we restrict it to 

I oo 



sequences {Nn}'^=i C N such that N'^M is an integer for each n G N. 
(3) If X, f, and /at are as in part 4.13| and a > (3 as in Theorem 



4.9 



then h% + AA^' 



h% + X\\--f\\L^j2) asN 
for /i^ + AiV-2| 



oo in the L-'^(T^) topology. A compactness result as in Theorem 



4.10 



J- I ^ 

Ml 

holds 



N\l- 



If instead the domain of h^ is restricted to for a fixed M G [0, 1], then the T -convergence and 



compactness results from Theorems 4.9 and 4.10 remain valid, with the domain of /ioo,o restricted to 
BVm{T^;{0,1}). 

We give a sketch of the proofs. 



(1) This follows directly from Theorem 3.6 



(2) For the fidelity term first we note that for v, f^ e An, 



N 



N\ 



Hence, ioic v £ An, u £ L^'(T^), / G C(T^) and fN £ An the discretization of / on Gat we have 



< \\v - fN 



< / \v-fN\- 
?^ + /llip(T2) < \\V 



/II 



Ilp{T2) 



+ 11/-/. 



N\ 



LP (T2)- 
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Since (by a Taylor series argument) /tv — >• / in LP(T^) as — t- oo the sequence of fidelity terms is 
continuously convergent and thus by Lemma |2.5| the F-limit follows. Clearly the compactness isn't 
harmed (even helped) by adding an extra term to the functional. 



For the mass constraint, since u G the conditions on M are necessary as explained in Section 



2.5 



L^(T ) convergence preserves average mass and hence the constraint is compatible with (LB) and 
compactness. For (UB) the recovery sequence {un}^^i used in the proof of Lemma 4.6 has support 
contained in the support of u and hence if Jj2 u = M, then Jr|j,2 Un < M. We can construct a similar 
recovery sequence {un}^=i, where Un has support on all grid squares which intersect suppn. This 
sequence satisfies all the required properties of a recovery sequence and has /.jj.2 Un > M. Hence, under 
the assumption on {Nn}^^i which assures that the mass condition can be satisfied for each Nn, there 
exists a recovery sequence {un}^=i where n„ takes the value 1 on suppn„ as well as on a select chosen 
number of squares which lie in suppu^ \ suppu„. For these combinations of {Nn}^^i and M (UB) is 
compatible with the mass constraint as well. 

(3) For the fidelity term we can use the same arguments as above. 

For the mass constraint we note that now Un G -^Nn and the limit function u G BV{T'^; {0,1}). 
This means each choice M G [0, 1] is allowed in the mass constraint. As above, because of the 
convergence this mass constraint is compatible with (LB) and compactness. For (UB) we note that the 
proof of Lemma 4.12 followed Lemma 4.6, so our argument here is very similar to that for h^fi above, 
with the added bonus that we do not need to restrict ourselves to specific combinations of {A''„}^^ 
and M. Let the recovery sequences {un}^=i and {un}^=i be as above. Now construct another recovery 
sequence {un}^=i by setting Un = Un on suppn„ and Un = CnUn on (suppu„)'^, where c„ G [0,1] is 
chosen such that for each n the average mass constraint is satisfied. 



r-limits for the the discretized 
Ginzburg-Landau functional k^^g 



1.5 



In this section we will study the convergence properties of k^^^ from 
N oo and then e — )• 0. The resulting F-limits are given in Section 
obtained by substituting e = N^" for well chosen a > and then considering the limit N 



5.1 



We first take the limit 
The simultaneous limit, 



in (2.3) is studied in Section 5.2 



oo for k'^ 



There has been a series of (recent) papers dealing with convergence of discrete energies to integral 
energies among which [T71 [5l [HI El [71 [8], all expanding on the ideas in [15, Chapter 4]. For many of 
the proofs we will use ideas from [5] in which the authors prove a F-limit of integral form exists for a 
general class of grid based functionals. Here we study a specific functional and hence can prove more 
explicit results. 



Remark 5.1. In Sections 



5.1 



and 



5.2 



we prove the results for , but they can easily be generalized to 
T'^ for any d G N, if we let the scaling factor in the first term of /c7v,e in (1-5) be eN'^~'^ instead of eN^ 
and the factor in the second term e~^N^'^ instead of s~^N^^ and we change kf, in (5.2) accordingly: 



N N 

i,j=l i,j=l 
N N 

k%{u) :=iv2-^-- ^(ui+ij - Uijf + {ui,j+i - u^^jf + N^-'' W{u^,j). 

i,j=l i,j=l 



(5.1) 
(5.2) 
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The admissible range of a in Theorems 5.6 and 5.7 is dependent on the dimension d in a way which 



will be made precise in Remark 5.10 For the extra fidelity term in Theorem 5.12 the scaling factor 



then needs to be instead of ^. 

5.1 Sequential F-convergence and compactness: first — )■ oo, tlien e — )■ 

We will prove kN,s —> koo,e as N ^ oo, where koo,e is defined for u G L^(T^) as 



''oo,e 



U 



e |Vn|2 + £-1 W{u) ifuG ^^^'^(T^), 
+0O otherwise. 



We see that /coo.e is the Ginzburg-Landau functional from (1.1). As explained in Section [l.l| it is known 



that this functional F-converges in either the L^(T^) or L^(T^ topolog}j^as e — )• to the total variation 



(1.2|). To be precise its F-limit is 

■oofiiu) :-- 



a{W)J^,\Vu\ iftxG W(T2;{0,1}), 
+00 otherwise, 



where a{W) := 2 yW{s) ds > is a constant depending on the specific form of VF, in particular on 
the transition between its wells, [15]. The sequence of functionals is equi-coercive as well. 

For fixed u G C^(T^) we have pointwise convergence (in C^(T^)) kN,e{u) — )• fcoo,e('u) as A'' — )• oo 
(more details below), but the dependence of the discretization errors on derivatives of u prevents us 
from concluding uniform convergence. F-convergence offers a useful middle ground between pointwise 
and uniform convergence and can thus be seen as an extension of classical numerical analysis results. 
The pointwise convergence follows from combining the boundedness of with the trapezoidal rule 
(for fixed i) 



T 



fix, y) dy = V f{x,j/N) + N-^-^{x, 0, for / G ^^(T^), 



for some /-dependent ^ G (0, 1), and the finite difference approximation of the derivative 

— {i/N,j/N) = N[u{{i + l)/A^,i/iV) - u{i/N,j/N)\ (5.3) 
-\N-'^^,{{^ + n)/N,J/N) 

for some G [0, 1]. We see the dependence of the errors on derivatives of the functions / and u. 
We prove the following F-convergence and compactness results. 

Theorem 5.2 (F-convergence). k]\j^i, — )• k^o^e as N ^ oo in the L1(T2) orL^{T^) topology. 

Theorem 5.3 (Compactness). Assume W satisfies (VF2). Let {Nn}'^=i C N satisfy Nn 00 as 
n — 7- 00 and let C //^(T^) be a sequence for which there is a constant C > such that for all 

n G N kiy^^i;{un) < C. Then there exists a subsequence {iin'}5v=i {^"IJ^Li ^^^^^ a u £ W^''^{T'^) such 
that Un' u in L'^{T'^) as n' —)• 00. 



^Results are usually stated in the topology, but for example [51] \T5[ [25] note that the results can be stated in the 
topology as well. 
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The proof of F-convergence adapts the ideas that are developed in an abstract general framework 
in [5] to our situation. In Section 4.1 we constructed a bijection between Vn and An from (4.1). In 
what follows we will identify u G Vn with its counterpart u £ An and drop the tilde if no confusion 
arises. 

For u G An C -L^(T^) we define pointwise evaluation by identifying each u with its representative 
z G we define the difference quotients as 



which is piecewise constant on the squares S^j^ and for which pointwise evaluation is well-defined. For 



D'}^u{z):=N[u{z + ek/N)-u{z)], k £ {1,2}, (5.4) 

where denotes the k^^ standard basis vector of M?. 

With the identification between Vat and An we extend the functional to all u G L^(T^) as follows 

kj^ (u) = l ^ [{Dh^f + iD%u)^] + £-1 f^, W{u) if n G An, 
'"^ 1^ +00 otherwise. 

In the proof of the liminf inequality we will use the slicing method, [15\ Chapter 15], [16], which uses 
the following notation. Remembering ^ [0, 1)^ we define Jj := {(0, y) G : 3t G [0, 1) : (t, y) G T^} 
and Tl := {(x,0) G : 3t G [0,1) : {x,t) G T^} and for (0, y) G Tf, (x,0) G T| we define the sets 
Tj^y := {t G [0, 1) : (t, y) G T^} and := {t G [0, 1) : (x, t) G T^} and the functions ui^y{t) := u{t, y) 
and U2,xit) := u{x,t) on T^^y and T^ .^. respectively. 
In what follows e > is fixed. 

Lemma 5.4 (Lower bound). Let u G ^^(T^) and let {un}n=i C ^^(T^) and {Nn}n=i C N be such that 
Un u in L^(T'^) and Nn — )> oo as n —)• oo. Then 

koo,s{u) < liminf fcAr„,e(nn). 
n— ^-oo 

Proof. This proof is an application of arguments in [5, Proposition 3.4]. 

First consider the case where u G W^''^{T^), then we can assume without loss of generality that 
{^N„,e{un)}^=i is Uniformly bounded and thus u„ G An„- For {x,y) G S]^^ we define 

vi{x,y) := Un{i/Nn,j/Nn) + D\j^Un{i/Nn,i/Nn) ■ {X - i/Nn), 

vl{x,y) := Un{i/N„,j/Nn) + D%^Un{i/NnJ /N^) ■ {y - j/Nn). 

In what follows k G {1,2}. Note G BVi^"^). We denote the densities of the absolutely contin- 

and 

N 



uous (with respect to the Lebesgue measure) part of the measures Vxv}^ and V^u^ by and 
respectively. Then for {x,y) G {S]^)° we have 



-^(x,y) = D}^^un{i/Nn,j/Nn) and -Q^{x,y) = D%^Un{i/Nn, j /Nn). 
We deduce — t- u in L"'^(T^) as n — t- oo from 
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The latter term converges to zero by assumption. The former we bound by \\v^ — UnWi'^ij'^) using 
Holder's inequality. We then note that from the uniform bound on {feAr„^e(«n)}^i we have 



^" r i2 
Y,[D%Un{ilNn,j/Nr,)\ <CNl 



for some C > and hence (if fc = 1; similarly for k = 2) 



[ (vl - Unf = V /" [Dlf^un{i/Nn,j/Nn) (x - i/N^)] ^ dy dx 

/t2 ^ — ' /q''J 

^" /■(i+l)/Afn (-i 

= J2 [Dh,Mi/Nn,j/Nn)] / {X - i/Nn? dx < -N'^ 

i,j = l -^V^n ^ 

For H^-a.e. y G the slice G W^''^{T1 ). By Fubini's theorem and Fatou's lemma 



lim inf / ( ^^—^ \ = lim inf 



V dx 



\{vl)'fdtdy 



> I lim inf / \{vl,)[ J"^ dt dy , 



Because 



lim inf / — ^ < lim inf kN„,e{uNj < oc 

n— s>oo Jj2 \ OX / n— !>oo 



(5.6) 



we have that, after possibly going to a subsequence, for T-L^-a.e. y € the sequence {||('t^n)'i,j/llL2(T2 )}-^i 
is bounded and hence weakly convergent in L^(Tf ,^). Since for Ti^-a.e. y G we have (f^)i,y — )■ ui^y 
in L^(Tf j^) we identify the limit as {v^)[^y u\ y in L^(Tf ,^), for T-L^-a.e. y G Tf (see Lemma 



B.l 



m 



Appendix |B| for details). Hence, by the weak lower semicontinuity of the L norm, Fatou's lemma, and 



the completely analogous results for U2^x and we get 



dy 



lim inf 



> 



\ui dt dy + 



\ dx ' \ dy ' 

Putting the slices back together using [32, 4.9.2 Theorem 2]) we deduce that 



|ti2 dt dy. 



lim inf 



\ dx 



\ dy J 



> 



(5.7) 



For the other term in the functional we use Fatou's lemma, the continuity of W and the a.e. pointwise 
convergence of Un to u to find 

liminf [ W{un)> I liminf H^(n„) = f W{u). 

n— s>oo Jr^2 Jj2 n-^oo Jj2 

This proves the result for u € W^'^{T^). 
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Now consider the case where u G L^(T^) \ l^-'^'^(T^). Assume that 

n—>-oo 

then after possibly going to a subsequence for each n G N G An„ and the sequences { Jj2 iD^^Un)'^^'^_^ 
are bounded. We can follow the same slicing method as applied above up to equation ( |5.6[ ). Again 
we find that for T-L^-a.e. y G Tf the sequence {||(^^n)i,i/llL2(x2 )}'^=i is bounded and combined with 
{vDi^y ui^y in L^i^ly) for n^-a.e. y G Tf we deduce that m^y G W^''^(Jly) for Ti^-a.e. y G Tf 
(see Lemma B.l in Appendix |B] for details). We then continue as above to arrive at (5.7) and conclude 
that u G W^'^Jt'^), which contradicts the assumption that u VF^'^(T^). Hence liminf /cat^ eiun) = oo, 

n— >oo ' 

which concludes the proof. □ 



oo as n 



oo and 



Next we prove (UB') (see Section 2.4). 

Lemma 5.5 (Upper bound). Let u G L^iJ"^) and {iVnj^i C N 6e such that iV„ 
let fc" be the T -upper limit ofk^^^e as n ^ oo with respect to the L^(T^) topology, then k'^{u) < koo,s{u). 

Proof. This proof is an adaptation of the ideas in [Sj Proposition 3.5]. 

The case where u G -L^(T^) \ ly^'^(T^) is trivial. For the case where u G PF^'^(T^) we first assume 
that u G C°°(T^). Define a sequence {un}'^^i such that for each n G N tin G -4iv„ in the following way. 
If {x,y) G S]^^ then n„(x) := u{i/Nn,j/Nn). Then (by a Taylor series argument) ?i„ — )• ti in L^(T^) as 
n oo. Let z/]^^ := {i/Nn,j/N„) G Gat^, then 



Jensen's inequality then gives 

(Z?jv„^n(i^^l))' = 



1 a-u 



9m . 
dx 



{v^ + sei/Nn) ds 



Qj^Ni+^^^/Nn)ds] < 



ds. 



Because u is smooth the Taylor series with remainder gives for z G 



du 
dx 



du 



du 



+ sei/Nn) = qz(^ + sei/Nn) + V— ((1 - c)x + cu'^J ■ ( 



dx 



for some c G [0, 1]. Using the fact that u and all its derivatives are bounded we find for some constants 
Cu > and Cu depending only on u 



-2 



^ / du ■ ■ 



ds 



( 



du . 



\dx 



{u^ +sei/Nn)] dsdx 



< 



< 



1 r /du 
^\dx 



(^Yiz) dz ds + Cu {Nn)-' + Cu {Nn 
Analogous estimates hold for D'j^^Un{y^^^)- 



z\+CuW'' 



ds dz 



r-CONVERGENCE OF GRAPH GiNZBURG-LANDAU FUNCTIONALS 



27 



Note that the sets S*^^ + sek/Nn, k G {1,2}, are just the squares shifted a distance s/Nn over 
the coordinate axes. Because we are working on the torus we get for some Cu > depending only on u 



< 



ds 



+ Cu{{Nnr^ + {Nn)-^)= [ |Vn|2 + C„((iV„)-i + (iV„)-2)_ 

We deduce that 

hmsup / {Dlf^Un) + {D%^Un) < / |Vnp. 

Because u G C°°(T^) n is bounded on and hence {|un|}^;^ is uniformly bounded on T^. Therefore 
there exists C > such that for all n G N W{un) < C. Because the constant C is integrable on T-^ 
we can use the dominated convergence theorem (or the reverse Fatou's lemma) and the continuity of 
W to deduce 

limsup / W{un) < / limsupVF(w„) = / W{u) 

(or use Jj2 {W{un) — W{u)) dx = W{s) dsdx < C fj2 \un — u\ dx). 

Combining the two inequalities above leads to 



lim sup /CAr„,e(u„ 
n— >oo 



< £ lim sup / [Dl^^Unf + {D'li^Unf + e limsup 

n— ^oo Jf'^ n— >oo 



J2 



< e J jVu\^ + e-^ I Wiu) = k^,eiu). (5.8) 



In the terminology of F-upper limit of Section 2.4 we have proven that k'^ ^{u) < koo,e{^) for u G 
C°°(T^). Since C°°(T^) is dense in l^^'-^(T^) (using VF-^'-^(T^) convergence) we use the lower semicon- 
tinuity of the F-upper limit to conclude (UB') for all u G ^"^'^(T^) as follows. Let {un]n=i C C°°(T2) 
be a sequence such that u„ — )• n in I^^'^(T^) as n — )• oo, then it also converges in L^(T^), hence 

k'l{u) < liminf /c''(n„) < liminf koo e{un)- 

Up to taking a subsequence Un ^ u pointwise a.e., hence W{un) — )■ W{u) pointwise a.e. Thus by 
possibly redefining u on a set of measure zero for n large enough we have that W{un) < W{u) + C for 
some C > 0. We can assume koo,s{u) < oo, hence W{u) is integrable on T^. Now we can again use the 
dominated convergence theorem or the reverse Fatou lemma and the continuity of W to find 

limsup / W{un) < / limsup PF(u„) = / W{u). 
Since Un ^ u in VF-^'^(T^) we have 

limsup / |Vii„p = / |Vu|^, 
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hence, 



k'eiu) < hminf koa^e{un) < hm 

sup /Loo,e (Wn) ^ ^oo,e(^)* 



□ 



Proof of Theorem 5.2| , Combining Lemmas |5.4| and 5.5 proves the F-convergence result. Note in 
particular that we have proven the lower bound for sequences converging in L^(T^) and the recovery 
sequence for the upper bound converges in L^(T^), hence we can conclude F-convergence in both 
topologies. □ 

Using a technique from [31'j 5.8.2 Theorem 3] we also get compactness for k^^e- 
Proof of Theorem 5.3, In what follows k S {1,2}. By (5.5) we have for all n G N u„ G An^ and 



E [ [{Dj^^unf + {Dl^unf] + [ W{un) < C. 
By assumption (VF2) on W we find that {||itn||L2('ir2)}^^ is uniformly bounded, hence there is a sub- 



sequence of {un}^=i (again labelled by n) and a u £ L^(T^) such that m„ ^ u in L^(T^) as n — )• 00. 
Moreover we see that {||-D^^u„||j;^2(']r2)}^]^ is uniformly bounded and hence there is a further subse- 
quence {tin'}^=i C {un}'^^i and awe L^(T^;M^) such that D'^ ^ Wk in L^(T^) as n' — )- 00. Let 
(/> G C~(T2), then 



Un'D^ ^(j) = Nn' / Un'ix)[(p{x + Ck/Nn') - (j){x)]dx 
T2 " 712 

= Nn' / [Un'{x - ek/Nn') - Un' {x)](l){x)dx 



By (5.3) we have /^a {D% ^(I)-V(I)- Ck) = KN~i^, for some K > hence D% ,0 ^ V(/>-efc in ^^(T^) as 



n' ^ 00. Combining this strong convergence for the difference quotient of (p with the weak convergence 
for Un' and its difference quotient we deduce 



uVcf) ■ Ck 



lim 



T2 



J2 



lim 

n'— >oo ./-jpa 



(j)Wk. 



T2 



Hence |Vm| = G L'^ij'^). We conclude that u G ^"^'^(T^). 

Finally, the strong convergence Un' — )• m in L^(T^) follows from the bound on Un by a discrete 



version of the Rellich-Kondrachov compactness theorem (see Lemma B.2 in Appendix^for details). □ 



5.2 Simultaneous scaling F-limit for k 



N 



In Section 5.1 we studied the F-limits of k^^g, by first taking N ^ 00 and then e — )• 0. We will now 



show we can take both limits at once if we scale e correctly in terms of A^. This is particularly relevant 
for numerical applications. We set e = N~°' for some a G (0, where q is the degree of polynomial 
growth of W' in condition (W4) and take the limit — >• 00 of k'^ in (5.2). 

Note that in contrast to the case for /i^ the order of the limits e — >• and N ^ 00 are reversed and 
thus we have an upper bound on a instead of a lower bound. 

We prove a compactness and F-convergence result. 

Theorem 5.6 (Compactness). Assume W satisfies (W3) and (W4) for given p > 2 and q > and 
a G (0, ^^). Let {Nn}'^^i C N satisfy Nn -^00 as n ^ 00 and let {un}'^=i C L^(T^) be a sequence for 
which tfiere is a constant C > such that for all n gN k^ (un) < C. Then there exists a subsequence 



{un'}^=i C {un}n=i and a u G BV(T^, {0, 1}) such that 



u in L\T^ 



as n 



00. 
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Theorem 5.7 (F-convergence) . Let W satisfy (VF3) and {W4) for given p,q > and assume a G 
(0, Then k% ^ fcoo.o as N ^ 00 in either the L'^{T'^) or L'^{T'^) topology. 



Using the difference quotient notation from ( |5.4[ ) we can mimic (5.5) and write 

j^a . )^ f N~-j^,[{D\,uf + {Dluf]+mi^,W{u) -line An, 
^" 1^ +00 otherwise. 

The proofs of this section make repeated use of the Modica-Mortola results |471 HSj H6l [5^ which 

for F^^^ ^ F^^ as 00. Note that condition 



show {L} and L^) compactness and convergence^*^ 



{W^) with p > 2 on the double well potential W is needed for the compactness result to hold (see e.g. 
[5H Proposition 3]). 

Let {un}'^=i C -L^(T^) be a sequence such that fc^ (?x„) < C. Below we will 



5.6 



Proof of Theorem 

prove the claim that there is a sequence {vn}'^=i C PF^'^(T^) such that \\vn — Un\\L2{j2^ — )• as n — )• 00 
and 

k%„{un) > F^Livn) + Rn, hm Rn = 0. (5.9) 

Given the veracity of this claim, it follows by the Modica-Mortola compactness result for F^^^ |47l H5l 

HSl [5l! that there is a subsequence of {vn}^=i (again labeled by n) such that u„ — )• u in L^(T^) for a 
u G BV{T'^, {0, 1}) (here we need condition (VF3) on W with p > 2). Using the triangle inequality and 
Ibn — Un\\L'2(^j2-^ — )• as n — )• 00 we then conclude that there is a subsequence of {un}^=i converging in 

L2(T2) to U. 

To prove the claim, first we show that 1 1 Un 1 1 00 (x2) = ©(A^n^^)- Assume not and let 7 > ^, then 
there is a subsequence (labeled again by n) such that for each n there is a square S]^^ on which 
\{un)i,j\ > A^'J. For definiteness assume {un^ij — A'^. By the uniform bound on /c^ {u^) we have 

N-'^iN^ - {un\+i,,f < C, (5.10) 

hence Uj+ij- = Q{N2). By induction over all the squares S]^^^ we find that Un = Q{Nn). Therefore 
II^".|Ilp(t2) = Qi^Tp), but by the coercivity condition (W3) on W (for any p > 0) the uniform bound 
on kfj- (un) demands 



Cl\\Un\\lr,^j2) < I W{Un)<CN-^ 



J2 



which is a contradiction. 

1 T- n ill T~i /-\ (-* T* 1 r ^-v* t x /-v I 'I- t f-t • l-J /~\ "i • / /i i \ f W_ 



Now for any u„ let be its bilinear interpolation: For (x, y) G 5](r'' define 



Vn(x,y) := nI {un)ij (^-j;^ 



^°As remarked in an earlier footnote, results are usually stated in the topology, but for example |51l 1151 125| note 
that the F-convergence can be stated in the topology as well. Compactness in follows from compactness in 
combined with the binary nature of the limit function. Finally note that the original results on bounded domains are 
easily adapted for the torus. Compactness is not hindered by the periodicity because no reg ularity beyond ^^(T^; {0, 1}) 
is needed for the limit. The lower bound generalizes immediately by restriction to sequences of periodic functions. The 
important properties of the recovery sequence for the upper bound are local properties near the boundary of supp u and 
so are also satisfied on a periodic domain. 
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+ K),+i,(x-^)(^-±l-j/) 



+ - ^) (2^ - it) + - i^;;) - it 

where for notational convenience wc have identified Un G '^Nn with its counterpart in Vn„- Thus 
defined Vn is continuous and ||i'n||L°°(T2) = ||'f^n||L°°(T2)- A straightforward computation shows 



{{_Un)i,j+l - (n„)i+ij+i)' 



First we note that there is a C > such that 



{{Un)i,j - iUn)i+l,j+iy 



< C 



((^n)i+lj — {Un)i,j+iy 



< c 



{{Un)i+l,j - iUn)i,jf + {{Un)i,j - 



Next we use periodicity to deduce 

XI - {Un)i+l,j+l) = ^ {{Un)i,j - , 

i,j=l i,j=l 

and analogously for similar terms. Using the uniform bound on A;^^ we find 
f r 2 

/ (Vn - Unf < CN-^ X [{{Un)i+l,j - {Un)i,jf + {{Un)i,j+1 - 

Thus, \\vn — UnW^^iji) — > as w — > oo. Another computation gives 

/ (^) [{{Un)i+l,j - {Un)i,jY + (K)i+ij+i - {Un)i,j+lY 

+ {{Un)i+l,j - {Un)i,j) {iUn)i+l,j+l - {Un)i,j+l) 

Using periodicity as above in combination with the inequality 

[{Un)i+l,j - {Un)i,j) {{Un)i+l,j+l - {Un)i,j+l) 
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< ■^(^{{Un)i+l,j - {Un)i,jf' + {{Un)i+l,j+l - {Un)i,j+lf'^ 



we deduce 



and the analogous result for j,^2 (i^^)^- 
Finally we note that 



dy 



K / W{Un)=K / W{Vn)+Rn, 



where 



Rn = K / (W^K) - WiVn)) < KM^'WUn " ^^nllLlfT^) (5.11) 
./t2 



Here 



M^':=max max \W' [s)\. 

xeT2 se[?;„(x),-u„(x')] 



By construction, since ||^in||Loo(']r2) = 0{Nn^'^), we have 



\Vnh^^j2^=0{N:/^). 



Hence, the maximum over s is achieved for some s = 0{Nn^'^). By the regularity of W and its 
polynomial growth condition {W4) we then have M^' = 0{Nn'^^'^), hence i?„ — )• as n — )• 00 by the 
choice of a. 

The inequalities above prove the claim and hence finish the proof. □ 



Remark 5.8. Clearly, the L (T ) convergence in Theorem 5.6 can be replaced by L (T ) convergence 
if desired. 



Lemma 5.9 (Lower bound). Assume W satisfies (W3) and (W4) for given p > and q > and 

2 



a e (0, J^)- Let u £ L^{T'^) and let {un}n=i C ^^(T^) and {iVn}^=i C N be such that Un ^ u in 



L^(T'^) and Nn — )• 00 as n ^ 00, then 

koofiiu) < liminf /c^ (u„). 

n— S.00 

Proof. Without loss of generality we can assume that A;^ (ti„) is uniformly bounded. In the proof of 



Theorem 5.6 we established that then estimate (5.9) follows, where the Vn are the bilinear interpolations 



of Un that converge to u in L^{T'^). Using the L-convergence result of Modica and Mortola [47 [ \Ab \ I^B l 
, specifically their lower bound, we find 

liminf fc^Jun) > liminf (f^^„(w„) + > k^fiiu). □ 



Remark 5.10. As noted earlier in Remark 5.1 our results (and proofs) generalize to T'^ if the terms 
in A;^ are rescaled properly depending on the dimension d. In this case we carefully need to reexamine 
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the admissible range for a in Theorem 5.6 and Lemma |5.9| (and hence by extension Theorem 5.7). In 
particular (5.10) becomes 

A^r'""(A^n^ - (^n).+l,i)' < C 

a + d-2 

and hence we need to choose 7 > ""'"^"^ and deduce that ||ti„||j;^oo(']i'2) = 0{Nn ^ ). Generalizing 
( 5. lip and the discussion that follows then leads to the conclusion that the admissible range of a is 
a G (0, ^"^"(^^3^^^ ) ■ In particular q in condition (VF4) should be chosen larger than d — 2. 

Lemma 5.11 (Upper bound). Let a G (0, 1), u £ -^^(T^), and {Nn}'^=i C N be such that iV„ cx) 
as n —)• 00 and let k" be the T-upper limit of k'^^ as n —)• 00 with respect to the L^(T^) topology, then 
k"{u) < koo,oiu). 

Proof. The case where u G ^^(T^) \ BViJ"^; {0,1}) is trivial, so assume u E ^^(T^; {0, 1}). First 
assume that suppw has smooth boundary Ssuppit. 

By the classical Modica-Mortola results used before /coo,e F-converges in both the L^(T^) and L^(T^) 
topology to /coo,o- Let {vn}^=i be the recovery sequence for this convergence with e = N~"-, see e.g. 
051 Proposition 2], [H §7.2.1], then each Vn G W^^'^{T'^) is a Lipschitz continuous function, — )• u in 
L^(T^) as n — )■ 00, and limsup A;^ ^-a{vn) < koo,Q{u). We denote the Lipschitz constant of Vn by Kn 

n— >oo ' " 

and note that Kn = 0{N") as n — )• 00. 

We now follow a similar line of reasoning as in the proof of Lemma 5.5 , but need to be more careful 
to deal with the lesser regularity of Vn (only Lipschitz continuous instead of C°°). For each i,j G /at 
and each n G N, define 



u^j^ := argmin 



where we used the difference quotient notation from (|5.4l). Note that since vn^ is continuous and S]^f is 



compact, the minimum is attained. Note that it may be that v]^^ G \ 5*]^^. Now define a sequence 

{un}^=i such that for each n G N G An„ in the following way. If x G S]^ , then Un{x) := Vn{v^^ )■ 
Note that by construction 



2 



First we check that ii„ — n in L^(T^). We estimate ||n„ — 'a||2,2(Tr2) < \\un — t;™ 11^2(^2) + \\vn — u\\i2(j2y 
We know that f „ — )• n in L^(T'^), for the first term on the right we use the Lipschitz continuity of Vn 
to compute 



Vn\\L2(j2) 



^ ^« E l^^l - ^ ^KlNlN~^N-\ 



Here we have used that there are nodes in the grid, the area of each square is N~'^, and 
for X G we have W]^^ — x\ < \/^N~^. Since = 0{N^) as n — >• 00 and a < 1 we have 

lL2(ir2) 



\\'i2(j2) — ;> as n — ;> 00 and hence it„ — ;> it in L^(T^). 
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Since vjy^ : — > M is Lipschitz continuous, if we fix eitlier y G T or x G T so are VNn{-,y) : T ^ M 
and vn^{x, •) : T — >■ M. Therefore by Rademacher's theorem the partial derivatives of vn^ exist a.e. on 
horizontal and vertical lines, hence we have for a.e. x G 



f dv f dv 

Dlfn'"n{x) = J -Q^i^ + sei/Nn) ds, D%^Vn{x) = J -Qy{x + se2/Nn) ds. 

By Jensen's inequality 



-^{x + sei/Nn)ds) < (^^{x + sei/Nn)) ds, 
and similarly for Dj^ Vn{x). For the finite difference terms in A;^ we now find, by construction of 



/ [{D},^Un{x)f + {D%^Un{x)f] dx 



E 



[{Dj,^Un{x)f + iDj,^Unix)f] dx 
[{Dl,^Vn{x))^ + {D%^Vn{x)f] dx 



£ E 



\-\dx 



ds dx 



\dx 



{x)Ydx+ f 



dx 



ds 



where we have used Pubini's theorem and the fact that we are working on a torus, so the union of all 
sets of the form S]^^^ + se^/Nn (for either = 1 or /c = 2) is the same as the Tmion of all S]^^, i.e., T^. 
To deal with the double well potential term in k'^^ we first note that for each x G and each 

n G N Vn{x) G [0, 1], hence W' is bounded on intervals of the form [vn{x) , Vn{i^'j^ )] by some C > 0. 



Therefore, for x G S^j^^, 



W{Un{x)) - W{Vn{x)) = W{Vn{iy%i)) " W{Vn{x)) 

< C\Vn{l^%{) - Vnix)\ < CKn\l/£ - X\. 



W\s) ds 



We thus find 



/ W{Un){x)dx= [ W{Vnix))dx + CKn V / |i 

< f W{vn{x))dx + CKnNlN-^N-\ 

JT2 



'Nr.. 



x\ dx 
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Hence fj2W{un){x)dx < Jrjpa VF(f„(a;)) dx + i?„,, where i?„ = 0{N^ ^). Since a < 1, Rn — as 



n 



oo. 



Combining both terms in k'^^ we find k'^^{un) < j^-a{vn) + Rn- We already know that 
Umsup A:^ < k^fi{u), so we have proved that for u € i?y(T^; {0, 1}) with smooth 5suppu we 

rn-oo ' " 



have hm sup k 



< koofliu) or in terms of the F-upper hmit: k"{u) < koofl{u). 



Now let u G BV{T'^; {0, 1}), not necessarily with the smoothness condition on the boundary. Then 
by O Theorem 3.42] there is a sequence {n„}^^ C BV{T'^; {0, 1}) such that each supptt„ has smooth 
boundary and lim fcoo o('Un) = kooo{u)- Hence by lower semicontinuity of the F-upper limit k" we 
conclude 

k"{u) < liminf k"{un) < liminf A;oo,o(^n) = koofliu)- □ 



Proof of Theorem 15.71 Since we have used L^{T^) convergence in (LB) in Lemma 
conv 
gies. 



5.9 



and L2(t2) 

convergence in Lemma 5.11 for (UB') we can now conclude F-convergence in either of these two topolo- 

□ 



5.3 Discussion of the range of a 

The range of admissible a in the results in the previous section is not only of theoretical interest, but 
is also important for computations. In simulations choosing e of the right order is a hard problem. If 
£ is too small in gradient fiow simulations this leads to the phenomenon of 'pinning', where the initial 
condition gets pinned down into the wells of W without changing its geometry. On the other hand, 
an e which is too large leads to immediate diffusion of the initial condition and loss of some relevant 
features. Our results do not directly address the gradient flow, but are in the same spirit. 

In the proof of compactness and the lower bound above we have assumed a G (0, where q is the 
degree of polynomial growth of W' . There are some reasons to believe this restriction could possibly 
be relaxed to a € (0, ,^), but we have not found a proof for this statement. 

First note that the dependence on q in the range of a comes from the discrete nature of the problem. 
Fundamentally it can be traced back to the lack of a chain rule for discrete differentiation. Trying to 
copy the classical Modica-Mortola result, we can deflne Wn := /J^" \/W{s) ds and estimate 



/ \VWn\ = N^-^ \{Wn)i+l,j - {Wn)i,j \ + |( 



Wn)i,j+1 - {Wn)i,j\ 



«J=1 



(«n)i + l,j ^ 2 / r(Un)i,] + l \2\\ 

^W{s) dsj + ( / ^W{s)ds] 



Nn 
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where G [{un)ij, {un)i+i,j] and G [{un)ij, (u„)ij+i] come from the mean value theo- 

rem. If we had control over the behavior of W in between grid points, we could use Cauchy's inequality 
to bound the expression above by kf^^{un)- Without condition (VF4) this control is lacking. We could 
do without this condition if we would somehow have an a priori L°° bound for the sequence {unj^i- 
This reflects a similar situation in the continuum case |54t Remark 1.35] where condition {W3) with 
p >2 can be dropped from the assumptions needed for compactness if an a priori bound is avail- 
able. By construction a uniform bound on Hn^H^oo^Tra) gives a similar bound on ||?;„||^oo(-']p2), hence 
under such an a priori bound we could drop conditions (W3) and (W4) from our assumptions and 
eleminate q (i.e., q = 0) from the restriction on a. The difference with the continuum case is that in 
our case (in the absence of an L°° bound) we need control over W and its derivative W . The scale of 
the discretization, A^~^, should be fine enough to resolve the variations in W' . 

Next note that in the proof of the upper bound we only use a G (0, 1). This restriction has a natural 
interpretation: If we interpret as the discretization spacing, then e = N~'^ for < a < 1 tells us 
that the discretization should be fine enough to 'resolve the diffuse interface' which, in the continuum 
case, has width of order e. The following example shows that for a > 1 the lower bound fails. Let 
u G BV(T'^, {0) 1}) be equal to zero on half the torus (say on [0, 1/2] x [0, 1)) and equal to one on the 
other half and let {un}^=i be a sequence converging to u in L^(T^) obtained by simply discretizing u 
on the grid G]\f„ for a sequence {Nn}^^i, — )■ 00 as n — )■ 00. Then W{un) = for all n. The finite 
difference term in k'^ {un) only has nonzero contributions along the boundary between the parts of 
the torus where u = and n = 1. Thus there are 2N jumps of order 1 and hence A;^^(u„) = 2N}~'^. 
If a > 1 this converges to zero, but the limit functional koo,o{u) > 0, which contradicts (LB). 



Finally, note that in (5.11) our estimate is not sharp since we use Holder's inequality to go from 
\\un — VnWL^^ji^ to \\un — Vn\\L^(T2y If instead a bound \\un — Wn||Li(T2) = ^i^n^) could be proved, 
possibly using the uniform bound on kfj^ (un), then in (5.11 ) the condition on a relaxes to a G (0, ^^), 



which in the absence of q would reduce to a G (0,1). We therefore conjecture that a G (0, ,^) is 



in fact the natural restriction for a (on T^, see Remark 5.10 for a discussion about the range of a in 
general dimensions), or, if an a priori L°° bound is available, a G (0, 1). However, it might be the case 
that a bilinear interpolation is not the right interpolation to attain this bound. 

5.4 Constraints 

In this section we show that addition of a fidelity term or imposing a mass constraint are compatible 
with the three F-limits we established, i.e., N ^ 00 for A;jv^e, e — )• for k^o^e and N ^ 00 for k'^. 

The Modica-Mortola limit k^o^e — ^ ^00,0 as — >• 00 in the LP(T^) topology, p G {1, 2}, is known to 
be compatible with a mass constraint, e.g. \45\ Proposition 2], |54t Theorem 1], |15l Proposition 6.6]. 
Furthermore, since an LP(T^) fidelity term, p G {1,2}, is clearly continuous with respect to LP(T^) 
convergence, it is also compatible with the F-limit. The theorem below addresses the other two F-limits 

for kiy,e and kfj. 

Theorem 5.12 (Constraints). (1) kN,e + AA^~^| • —/nIp ^oo,e + /x2 I ■ —f\p for N ^ 00 in the 
L^lT"^) topology, where p G {1, 2}, A > 0, / G C^(T^) and fj\[ G An is the sampling of f on the grid Gn 
(f ! fN cLnd their norms can also be defined on subsets ofT"^ and Gj\[ as in Theorem 4.13, part 4.13). 
A compactness result for k^^e + AA^~^| • — /|p as in Theorem 



5.3 



holds. 

If instead, for fixed M G [0, 1], the domain of definition of kN,e is restricted to {i.e., V^^ from 



Theorem 3.6 on the grid G^), then the T -convergence and compactness results for — )■ remain valid. 
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with the domain of koo,o restricted to V 

I —JN\p —I 



M 



(2) k% + \N-'^\--fN\p feoo.o + A /j,2 \--f\l for N ^oo in the LP(T^ topology, where p £ {1,2}, 
A, /, and f]\j are as in part 



5.12 



-f\p as in Theorem 



and a G (0, 



and Remark 



5.6 



-?+3. 



5.8 



as in Theorem 
holds. 



5.7 



A compactness result for 



If instead, for fixed M € [0, 1], the domain of definition of k'^ is restricted to VfJ , then the T- 



convergence and compactness results for N 



remain valid, with the domain of /Coo,o restricted to 



We give a sketch of the proofs. 

(1) The compatibility of the fidelity term with the F-convergence and compactness follows as in 
the proof of Theorem 4.13, part 4.13 As in that theorem the mass constraint is preserved under 
LP{T'^) convergence and so is compatible with both (LB) and compactness^^ To show that the mass 
constraint is compatible with (UB') as well we need to check two conditions. First, the recovery 
sequence which was constructed (in the proof of Lemma 5.5) for u £ C^{T'^) should either satisfy or 
be able to be adapted to satisfy the mass constraint. Second, for each u G H^^'^(T^) there should be 
an approximating sequence {un}^=i C C°^{T'^) which has constant mass. The latter follows directly 
by the use of normalized mollifiers to construct the approximating sequence. For the former condition 
we follow an argument reminiscent of the proof that a mass constraint is compatible with the Modica- 
Mortola F-convergence result for the continuum Ginzburg-Landau functional, see e.g. [45t [Ml I15j . 
Assume that Jr|j,2 u = M for some M > 0. Because u is smooth it has bounded derivatives on T^, hence 



(using the notation S]^ from (2.2)) 



Un — U 



= f Hi/Nn,i/Nn) - U{x)] dx 

<Cuy^ \{i/Nn,j/Nn) - X\2 dx < CuN-\ 



for some constants Cm, Cu, depending only on u. Hence for each n G N there is a (5„ = 0{N^ ^) such 



that Un := Un + 6n satisfies /rj 



Y2 '^n 



M . For this new proposed recovery sequence we compute 



kN„,e{Un) = kN„,s{Un) + £ W [W{Un) " W{Un)] . 

By Taylor's theorem, for x G T^, 

W{Un{x)) - W{Un{x)) = VF'(c„(x))5„, 

where c„(x) G [un,Un + Sn]- Since u is continuous and hence bounded on T^, the sequence {un}^i is 



equibounded and hence, because W G C (M), W'{cn) is equibounded. Therefore, the sequence {u„}^^ 
is indeed a recovery sequence for (UB) with u G C°°(T^) and satisfies the mass constraint. 



(2) As in part 5.12 the addition of fidelity terms is compatible with the F-limit and compactness 



and the mass constraint is compatible with both (LB) and compactness. For compatibility with (UB') 



^Both the fidelity term and the mass constraint are not only compatible with the compactness result, but even help with 



concluding uniform boundedness of either I/^(T^) or L^(T^) norm and hence can replace assumption (W2) in Theorem 
when compactness with respect to the correct topology is considered 



5.3 
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again we check two things: First that the recovery sequence which was constructed (in the proof 
of Lemma 5.11) for u £ BV{T'^; {0,1}) with d suppu smooth either satisfies or can be adapted to 



satisfy the mass constraint and second that for each u G BV{T'^; {0,1}) an approximating sequence 
{un}^! C BV(T'^; {0, 1}) can be chosen for which dsupp Un is smooth and which has constant mass. 
The latter condition is satisfied if we use the approximating sequence as in |9i Theorem 3.42] and then 
introduce a small dilation, diminishing along the sequence, of the support of each Un so that the mass 
remains fixed (see e.g. \50\ Proposition 7.1]). The former condition follows in a similar way to the 
construction above in the case N ^ oo for kf^f ^. 



5.5 Gradient flow for /c^v.^ with constraints 

To minimize fc7v,e either under a mass constraint or with a fidelity term we can use a gradient fiow. 
First consider the latter case: kN,e,x '■= ^Af,£ + A| • — /ll- For u, f G Vn we compute gTad{kN^e,x){u) G Vn 
via 



d_ 
dt 



kN,e,xi'^ + 



t=0 



{grad{kN,ex)iu),v)\ 



and then set for all i,jGlN 
This leads to the equation 

dt 



du 



dt 



-{gTad{kN,e,x){u)).j. 



2e (uij - Uk,i) + e-^N-^W'im^j) + 2X{uij - fi/j 

■ {k,l)€Af{i,j) 



(5.12) 



where the set of indices of neighbors of (i, j) is given hy M{i,j) = {{i — l,j), {i + l,j), {i,j — l), (i, j + 1)}. 
The overall prefactor 4"*" comes from the factor d~j^ which is needed to cancel the factor j in the Vat 
inner product. Here we assume the weights in this case to be equal to 1 (on existing edges). Equation 



(5.12) is the discretized analogue of the continuum Allen-Cahn equation with data fidelity 

2eAu - e-^W'{u) - 2X{u - /), 



du 



which is the gradient flow of F^^{u) + A||m — /||^2(']I■2)• 
If instead of the fidelity term a mass constraint is imposed the term 2X{uij — fij) gets replaced by 
a Lagrange multiplier 

N 

To illustrate we show simulation results using a fidelity term with / the characteristic function of a 
square. We use a one-step forward in time finite difference scheme to discretize the time derivative, 
i.e., 



n+l . 
id 



u: 



A'^'dt 



2e E {ul^-uli) + e-^N-^W'{ul^) + 2X{ul^-h,) . 
{k,i)eAf{i,j) 

Here dt is the discrete time step and the superscript n labels the time step. We start with a random 
initial condition u^. We use the inner product structure on Vn corresponding to the unnormalized 
Laplacian (r = 0). Using the structure corresponding to the random walk Laplacian (r = 1) instead 
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Figure 1: (A)-(G) show snap shots of the gradient flow (5.12) using the parameters in the text, (h) 
shows the corresponding time evolution of A:7v,e,A('u)- 
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only gives an overall multiplicative factor | in the right hand side of the gradient flow in (5.12) and 
hence is effectively just a time rescaling leading to qualitatively the same behavior. 

In Figure [T] we use the following parameter values: = 100, e = 5, A = 0.1 and dt = 0.01. We use 
W{s) = s^(s — 1)^ for the potential and / is data prescribed to be 1 in a square region and outside 
that region. 



Note that W' satisfies the growth condition (W4) with q = 3, hence according to Theorems 5.6 



and 5.7 and e should satisfy the relation N~°' = e for an a G (0, 1/3). The combination = 100 
and e = 5 used in Figure [T] falls outside this range (a ^ —0.35), but the simulations still give a good 
result. Our theoretical results give a good guideline for choosing a (especially the upper bound), but 
in practice different values of a, and hence e, can produce good gradient flow simulations. We have 
chosen this larger value of e for our figures to show a stronger diffusion. 



6 The continuum limit of nonlocal means 



In this section we study the nonlocal means functional from (1.6) for a given fixed ^ E C°°(T ). 
Remember that the weights are ujl,n '■= s~ with a,L > constants (possibly depending on N) 

and dL,N defined in (1.7). We define the limit weights ujL,a,<^e,c £ -L°^(T^ x T^) as 



where £,c > and Se := {z G : \zi\ + \z2\ < £}. The limit functionals 5^ : ^^(T^) R and 
qL : L^T^) are 



9ooiu) ■■= / u}L,aix,y)\u{x) -u{y)\dxdy, 
aioiu) ■■= / i^e,c{^,y)\u{x) - u{y)\dxdy. 
We prove a F-convergence result. 

Theorem 6.1 (F-convergence). (1) If a,L > are fixed, then g^ ^ g^ as N ^ 00, in the L^(T^) 
topology. 

(2) If a = ^ for some c > and L is such that L/N £ as N ^ 00, for some £ G (0, 1/2), then 
9n gio as N ^ 00, in the L^(T^) topology. 



As explained in Remark 6.5 we do not have compactness in this case. 



Note that ^at is the functional /o from Theorem 3.1 where the graph G is the grid Gn and the 
choices x = N~'^ and a; = ujl,n have been made. Two main differences between this functional and the 
previous functionals on the grid Gn we considered are that the graph is now completely connected and 
the weights are not uniform over the edges. For the latter reason it is useful to introduce a notation 
for the space of graph weights on Gat. Given a weight function oj and nodes nij,nk^i G Vn we write 
^i,j,k,i ■= ^{ni,j,nk,i). Define 



Wn ■■= {uj -.Vn xVn ^ [0, 00) : for ah i,j, k, I G 1,^2 ujij^i = (^k,i,i,j]- 
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Completely analogous to the identification between Vn and An which was introduced in Section 4.1 
we can identify Wn with 

:= : X ^ [0, oo) : for all x, y G u{x, y) = u{y, x)}. 



Identifying ujl,n G VVat with the corresponding ujl,n ^ and u G with the corresponding u ^ A'' 



N 



we can write 



9N[u) 



Ij2jj2i^L,N{x,y)\u{x)-u{y)\ 



dx dy if u G 



+00 



ifuGLi(T2)\^^. 



We prove Theorem 6.1 in two steps. First we show that uniform convergence of the weights suffices 



for F-convergence of and then we show that uniform convergence. 

Lemma 6.2. Let {ijJN}'^=i such that ujn ^ Oind oj^ ^ oj uniformly as N ^ oo for some 
u) G L°°(T2 X T^). ForN(^n define the functional Qoo ■ L^{T^) -^R by 

goo{u) := I / uj{x,y)\u{x) -u{y)\dxdy, 



then gjsf — )• ^oo as N ^ oo in the L^(T^) topology. 

Proof. Let {iVn}^=i C N such that iV„ ^ oo as n ^ oo, u G L^{T'^) and {n„}^=i C L'^{T'^) such that 
Un ^ u in. L-^(T^) as n — )■ oo. Assume that for each n G N n„ G A^j^^. Then 



^N„{x,y)\un{x)-Un{y)\dxdy- I I u){x,y)\u{x) - u{y)\dxdy 
{^N.a{x,y) - uj{x,y))\un{x) - Un{y)\dxdy 



T2 JT2 

+ / / ^{x,y){\un{x)-Un{y)\-\u{x)-u{y)\) dxdy 
<Ii+ / u}{x,y)\un{x) -Un{y) -u{x) + u{y)\dxdy < I1 + 2I2, 

Jt2 Jt2 



where for simplicity we have used the notation 



h{u) := I I \ujNr,{x,y)-u:{x,y)\\un{x)-Un{y)\dxdy 

/t2 712 



hiu) 



T2 JT 



oj{x,y)\un{x) — u{x)\ dxdy. 



Since cjjv^ — )• a; uniformly as n — )• 00, there is a sequence of constants C„ > such that for n large 
enough |a;Ar,^(x,y) — uj{x,y)\ < Cn and Cn —5- as n —)• 00. Hence, 

h<Cn / - Un{y)\ dxdy. 

Jt2 yT2 

Because u„ G we have \un{x) — Un{y)\ < 2 for almost all {x,y) G x and hence /i — 
as n — 00. Furthermore I2 < ||a;||j;^oo(']i'2) ||un — n||j;^i(']p2) and thus /2 — ^ as n — 00. We conclude 
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that Urn gN„{un) 

n— >oo 



u in L1(T2) as 



goo{u) for any sequence {un\^=i such that u„ G ^^^^ and Un 

n — >• oo. In particular (LB) is proven. To prove (UB') all that remains is to show that there exists such 
a sequence. 

First assume u G C°°(T^) and define Un{x) := u{i/ Nn, j / Nn) where i,j G In are such that x G S]^^ 
from (2.2). This sequence satisfies the required conditions, hence (UB') is proved for u G C°°(T^). We 



conclude the argument by using the lower semicontinuity of the upper F- limit and density of C°°(T^) 
in L^(T^) as in the proof of Theorems 5.5 and 5.11 to deduce (UB') for u G L^(T^). □ 



Remark 6.3. Note that qn does not converge uniformly to goo, because if li G L1(T2) is such that 
for all G N u A^j^, then \gN{u) — goo{u)\ = oo. If we would restrict the domains of gN and goo to 
continuous u and define g'Ar to be 



gN{u) := 



T2 JT2 



ujN{x,y)\uNix) - UNiy)\ dxdy, 



where un(x) 



u{i/N,j/N) where i and j are such that x G S^'-' 



then gN ^ goo uniformly by 



the estimates in the proof of Lemma 6.2 By [28, Proposition 5.2] gN then F-converges to the lower 
semicontinuous envelope of goo- 



Lemma 6.4. (1) If a, L > are fixed, then ujl,n '^L,o- uniformly as N ^ oo. 

(2) If a = ^ for some c > and L is such that L/N i as N oo, for some £ G (0, 1/2), then 
^L,N ^l,c uniformly as N ^ oo. 

We defer the relatively straightforward proof to Appendix [B| 
Proof of Theorem 6.1[ Combining Lemmas |6.2| and |6.4| the result follows directly. □ 



Remark 6.5. It is important to note that for gN we do not have a compactness result in the L^(T^) 
topology. If we have sequences {Nn}'^=i and {un}'^=i such that Nn — )• oo as n — oo and u„ G "4.^^, 

the bound on ^00(^2) allows us to conclude that Un ^ u for some subsequence (labelled again by 
n) and some u G L°°{T ). In order to deduce L^(T^) convergence we would need some information 
on the derivatives (or finite differences) which we do not have when the weights w are nonsingular. A 
uniform bound on gN„{un) adds no useful information since gNni^n) < 1 per definition, if m„ G 

A simple counterexample is the case where $ is constant, hence the graph weight function u = I. 
Let MTV £ be a checker board pattern on Gn, i-e., as function in {un)o,o = and {uN)i,j / 
{uN)i+i,j = {uN)i,j+i for ah ij, then 



gN{uN)= / \un{x) -UN{y)\dxdy = 2\suj)puN\\{sup-puNy\ <2. 
Jt2 Jt^ 

However no subsequence of {un}n converges in L^(T^) as can be seen as follows. Let M ^ N then the 
square S'J contains o((f )^) squares of size M ^ by M ^. On approximately half (at least 0(1)) of 

these UM + UN, so \un - um\ = NHr^o(^{^f)o{l) = 0(1). 
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7 Discussion and open questions 



In this paper we have shown various F-convergence results. The convergence of /e in Section [3] shows 
that we can extend the classical Modica-Mortola F-convergence result for the Ginzburg-Landau func- 
tional to graphs, if we are careful about the precise scaling. The discrete nature of a graph forces us 
to not include an e in the finite difference term, unlike the e in the gradient term for the continuum 
Ginzburg-Landau functional. As has been shown on a specific regular square grid in Section [4j this 
has consequences for the limit functional, which now behaves like an anisotropic instead of isotropic 
total variation. We recovered the isotropic total variation for the regular grid case in Section [5] by 
taking an approach reminiscent of classical numerical analytic results, instead of graph based results. 
Specifically, to do this we need to choose a scaling in line with standard finite difference and quadrature 
methods and make the limit N ^ oo dominant over the limit e — >• such that, in a sense, we first get 
back to the continuum case, before passing to the total variation. The lesson in here is twofold. On the 
one hand it shows that one has to be careful when discretizing on a grid not to pick up grid direction, 
which has been known to numerical analysts for a long time. On the other hand however, it entices 
us to look at graph based functionals and nonlinear partial differential equations in their own right, 
because they can behave in surprising ways if the topology of the graph is allowed to interact with the 
functional or PDE. This conclusion is reminiscent of the behavior which is found in [3T]. In that paper 
the authors study the limit of the graph Laplacian on a graph which is constructed by sampling points 
from a manifold. They find the limit is independent of the sampling distribution only for a specific 
scaling of the graph Laplacian. 

In Section [6] we studied the limit of a functional of nonlinear means type on graphs, showing that 
while the limit exists, the nonlocal nature of the functional leads to loss of compactness. This is not 
expected to be a specific problem of the graph based nature of the functional, but of the nonlocality 
and as such is expected to be present for a continuum version of qn as well. 



One question raised in Section 5.3, is whether the range of a under which F-convergence and 
compactness of kfj can be proven, can be extended to (0, 1). This is an important question in practice 
when running gradient flow simulations. A choice of e which is too small or large can lead to either 
pinning or too fast diffusion respectively. 

The F-convergence results for /i7v,e and /cat^e naturally lead to the question of the limit behavior for 
other graphs. In order to have a good interpretation for that question it is in the first place necessary to 
have a structured way in which to increase m, the number of nodes for the graph. A triangulation might 
be the natural next step. A random Erdos-Renyi graph [30j also carries a natural rule how to connect 
new nodes to the graph and may be an interesting exploration into the question whether the Ginzburg- 
Landau functional on a graph without explicit spatial embedding can possibly have 'continuum' limit. 
For arbitrary graphs it is less clear how to add new nodes in a structured way. One option could be to 
construct a sequence of graphs where in each next step each existing edge is bisected by a new node. 

A question that is very relevant for the applications of the Ginzburg-Landau functional is that of 
stability of minimizers with respect to perturbations of the graph {e.g. perturb the weights or add or 
delete nodes). For example, in data analysis, if the nodes represent data points and the edge weights 
measure similarity, it is quite likely that noise is present in the weights. 
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A Choice of Hilbert space and difference structure 

In this section we will give some background information and derivations concerning the choices made 
in Section 2.2 that defined our graph operators and functionals. 

We start by associating V and E with the finite dimensional vector spaces and M™('"~i)/2 
respectively. We will turn these vector spaces into Hilbert spaces by defining inner products on them. 
We follow the procedure described in [HI Section 2 For u, t> S V and ip, (p G £ we define 

{u,v)v := ^UiVia{di), ((^,0)^:=^ ^ (fij4>ij(3{ujij), 

where a, (3 : [0, oo) — )• [0, oo) are functions yet to be determined. Note that a priori we allow a and f3 to 
take the value zero, which means that the above 'inner products' might be positive semi-definite and 



not positive definite. We will get back to this issue after we have decided on our choices of a and f:^ 
As in [35] we also define the dot product for tp, (j) £ £ as 



• 4>)i -=7^^ ipij(l)ijP{u 



2 ^ r.JT.j:A^i. 

Note that (/9 • G V. As explained in ^41] we can now define the difference operator or gradient 
V : V ^ £: as 

(Vn)ij := 7(a;ij)(uj - Uj), 

where 7 : [0, 00] —J- [0, 00] is a third yet to be determined function. With this choice for the gradient 
we find that its adjoint, the divergence div : i5 — t- V, is given by |4H Lemma 3] 

This expression follows from the defining property of the adjoint: (Vu, ip)s = {u, div(/3)v for all G V 
and all 93 G <5. 

Now that we have inner products, a gradient operator, and a divergence operator, we can define 
the following objects: 

• Inner product norms ||u||v := \/ {u, u)v and \\<f\\e := \/ {(p, tf)s. 



^^We slightly deviate from [41] ■ Instead of sums YlTLi they use averages — JZ^^i, which is a choice not unanimously 
adopted in the literature, but leads to cleaner convergence statements in [JT]. We could adopt this convention in this 
paper, but it would not significantly alter our results, mutatis mutandis. 

^^Note that if a and /3 are such that positive definiteness is satisfied, these inner products do indeed turn V and £ into 
Hilbert spaces, since convergence with respect to the induced £ norm preservers skew-symmetry. 
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• Maximum norms -"^^ Hiillv.oo := max{|tij| : i G Im} and 
llv'lk.oo := max{|(^ij| : i,j £ Im}- 

• The norm corresponding to the dot product \\(p\\i := {(f ■ Lp)i. 
Note that || • ||£:,dot G V. 

• The Dirichlet energy ^llVuHl. 

• The graph Laplacian A := div o V : V — )> V. So 

(Au)i := — ^ V /3(a;ij)7^(wij)(ni - Uj). 

• The isotropic and anisotropic total variation TV : V — )• M and TVa : V — )■ M respectively: 

TV(n) := max{(div(/7, -u)-); : £ S, max H^^Hi < 1}. 

TVa(n) := max{(divv3,u)v : G <S, ||(/?||£-,oo < !}• 

We note that by the property of the adjoint we can also use (Vu, ip)^ in the definitions above instead 
of (divyj, u)v- In this finite dimensional setting these maxima over unit balls will be achieved, hence 
we are justified in using max instead of sup. 

Before we start making specific choices for a, /3, and 7 it is interesting to make some general 
observations which do not depend on these choices. 

• We can alternatively derive the Laplacian via the variational principle from the Dirichlet energy 
as follows. Consider u,v and t G M, then 



2 



\ ^ (3{u}ij)j'^{uJij){ui-Uj){vi-Vj) 



= ^ f3{ujij)^'^{ujij){ui - Uj)vi 

= 2^ X ' [Ui-Uj)via{di) = {Au,v)v. 

If we choose v = u this also shows that an analogue of 'integration by parts' holds: 

(An,n)v = ^ ^ I3{u}ij)-f^{u}ij){ui-ujf = \\\7u\\l. 

• We can also give a variational TV-type formulation of the Dirichlet energy itself via: 

||Vti||£: = max{(div(^, n)v : (p £ £, \\p>\\s ^ !}• 



^"^To justify these definitions and convince ourselves that there should be no /? or 7 included in the maximum norms 
we define ||<p||£p := IX^ijg/ VijPi^^ij)- Adapting the proofs in the continuum case in e.g. [TJ Theorems 2.3 and 
2.8] to the graph situation we can prove a Holder inequality ||</7<j!>||£,i < ||i/'||£,p||<^||£,q for 1 < p, g < 00 such that 

f 1 \ ^ ' 

i + i = 1, an embedding theorem of the form ||(^||£,p — ( 2 Pi^ii) | llvlk,? for ^ ^ P ^ 1 ^ 00 and the limit 



lim ||<,o||£,p ~ W^pWs.oo. a similar result holds for the norms on V. 

p — ^00 



r-CONVERGENCE OF GRAPH GiNZBURG-LANDAU FUNCTIONALS 



45 



To see why this holds we first remember that (Vn, ip)^ = (div(/?, n)v- Then we see that by the Cauchy- 
Schwarz inequaUty 

{Vu,ip)£ < WVuWsWipWs < llVullf. 

Equahty is achieved when 



ip = ^ (u) :-- 



||VU||£ II 11'^ ' ' 

if||Vu||f = ' 



which is permissible since < 1. 

We will now make particular choices for a, /3, and 7. Our choices will be driven by the desire to 
satisfy the following properties. 

(1) We will consider a family of graph Laplacians indexed by a parameter r < 1 (not to be confused 
with the p-Laplacian from the literature). As it turns out only the choice of a is influenced by the 
choice of r. The Laplacians we consider are 

(A,u)i := dl-'ui - ^ "^Uj =Y'^iui- Uj) 



As explained in Section 2.2 by choosing either r = or r = 1 we recover the unnormalized or 
random walk Laplacian respectively. To construct the symmetric normalized Laplacian as it appears 
in the literature requires a gradient of the form 

Ui 



{Vu)ij = 7(wij)(- 



/Ji 



{cf. |H]). This falls outside our current framework and hence we will not consider it here,^^ Some 
discussion of the pros and cons of different graph Laplacians can be found in e.g. I58j. 

(2) The Dirichlet energy is given by 2l|Vn||| = \ Yli j^i^n^^i^^^ independently of the choice 
of r in the Laplacian. 

(3) The isotropic total variation is 



T\{u) = l|Vn|k = \^ JYI ^ijiui-Uj)^ 

{cf. [35j where TV is called nonlocal TV because the graph is assumed to be embedded in an Euclidean 
space and so what is local on the graph (neighboring vertices) might not be local in the embedding 
space) . 

(4) We will consider a family of anisotropic total variations parametrized by the parameter q G 

[1/2, ir 



16 



2 ^ 



TV,g{u) = {Vu,sgniVu))£ = l Y 

The parameter q comes in via the definitions of /? and 7 and the signum function is understood to act 
element-wise on the elements of Vu. 

Let us consider all the points above to find out what conditions we have to impose on a, /3, and 7 
to satisfy this list of requirements. 



^^As a word of caution we note that the use of the symmetric normalized Laplacian in combination with a double well 
potential W with wells that are not symmetrically placed around (as in the case where the wells are at and 1) causes 
an asymmetry between the two phases that is typically unwanted. 

^®We can take g G E if we interpret oj^^ as zero whenever Uij — 0. 



46 



Yves van Gennip, Andrea L. Bertozzi 



1. As can be seen in the definition of tlie Laplacian (c/. also [41^ Definition 7]) in order to get the 
desired Laplacians we have to choose a, (3, and 7 such that for each LOtj and each df. 

a{di) d\ ' 

Specifically j3{ujij)^'^{ujij) = coij for any choice of r and a{di) = d\. We will see below that the 
choice of a is irrelevant for the points [2}|4] and hence all choices of r are compatible with what 
follows. 

2. For the Dirichlet energy we compute 

Since the graph Laplacian appears as the natural operator in the Euler-Lagrange equation asso- 
ciated with the Dirichlet energy it is not surprising that we do not get any extra conditions on a, 
/3 or 7 from the Dirichlet energy which we didn't already get from the Laplacian. It is interesting 
to note though that the Dirichlet energy does not depend on the choice of a (and hence r in the 
Laplacian) at all. 

3. For the isotropic total variation TV we use the Cauchy-Schwarz inequality on the dot product 
norm to get {Vu,(f)£ = X]jg/^(V?i • (f)i < Yyiei,n - ^ieim ll^^ll*- '^^ achieve equal- 

j£ II Vnll ■ 7^ 

11^^*11! ' . Again we do not require extra conditions on 

if llVnII,- = 



it} 
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/3 and 7. They will be determined by the last requirement 



18 



4. To compute TVa we use the bound on the maximum norm of (p to find 

Wi3\\Ui-U0{u)ij)-i{uJij) 

\ui-Uj\l3{L0ij)'y{L0ij). 

To achieve equality we can choose (p = ip"" := sgn(Vn), i.e., (pfj = sgn{uj — Ui) for i, j such that 
"f{oJij) > and (^"^ = otherwis^^ Hence 

TVa(n) = ^ J2 \ui-Uj\(3{ujij)j{ujij). 



^'^Note that demanding ip'^^ to achieve equality does not determine it uniquely on the set of vertices for which || VttHi — 0. 
^*If we would have defined the dot product {(p ■ (fi)i :— | X^jgj 'fiij4>ijd{iiJij) for a function 5 possibly different than 

/3, the requirement on TV would have led to the condition ^ = ^ij- Together with /3(cji_,)7^(aj,;j) — ojij from 

point [T] this gives 5 = /? as we have assumed all along. 

^^Note that we can change 1^9" on the set of vertices for which \7u = without losing equality. See also footnote 



17 
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If we now choose /3{ujij) = ufj ^ and jiiOij) = oj\- then TVa = TVag while /3 and 7 satisfy the 
necessary condition I3{ujij)'y'^{u)ij) = coij from the previous points. 

These choices for a, (3, and 7 lead to the inner products (or semi-definite sesquilinear forms), operators 



and functions presented in Section 2.2 



It is interesting to consider the conditions under which ||Vn||£- = and ||Vn||j = 0. 

\\Vu\\£ = 0<^ ^ ujijiui - Uj)"^ <^ V(i, j) G [uij = 0V Ui = uj]. 

This means that ||Vu||£- = iff li is constant on connected components of the graph. Similarly 



|Vu||j = 44> ^ ujijiui - Ujf 4^ Mj G Im [wjj = V lij 



u 



hence ||Vit||j = iff n is constant on the set {rij €z V : j = i \/ Cij G E} consisting of neighboring 
vertices of the i^^ vertex plus the i**^ vertex itself. 

We see that these conditions do what we would hope and expect them to do, even if the choice of q 
has made the f-sesquilinear form semi-definite, i.e., ||Vn||£- = gives global (per connected component) 
constants and ||Vn||j = gives local constancy. 



B Deferred proofs 



The next lemma was used in the proof of Lemma 5.4 



Lemma B.l. Let {un}^^i C L^{T) be a such that Un ^ u in L^{T) for a u £ L^(T^) and n'„ —^vin 
L2(T) for ave L'^{T). Then v = u' (and thus u G W^''^{T) ). 

Note in the proof below that this result in fact does not depend on the dimension and holds on T". 

Proof of Lemma EUl Let ip G C^{T), then 



= lim / Lp{v — u'jf) = lim / \ipv — ip'u — ipu'^ + (p'u] . 
There is a C > such that 

lim / \—(f'u — (pu'n] = lim / (f'[un — u]< lim C\\un — u\\li(j) = 0, 



hence we conclude j^^ipv = — ip'u. □ 

The next lemma is a discrete Rellich-Kondrachov type compactness result used in proof of Theo- 
rem [531 

Lemma B. 2 (Discrete Rellich-Kondrachov compactness result). Let {un}^i C L^(T^) and {Nn}'^^i C 
(0, 00) be sequences such that as n ^ 00 we have Nn — >• 00, ^ u in L^(T^) for some u G L^(T^), and 
the difference quotients {see (5.4)) \\D^ UnWi'^ij'i) {k G {1,2}) are uniformly bounded. Then Un ^ u 
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Proof. For e > let := JeUn € C°^(T^) be a mollified function on the torus, defined to be the 
solution to the heat equation after time e with initial condition u„: 

where u^{k) = Un{x)e~^'^^^'^ dx. We proceed in two steps. First we need to prove some properties 
of the mollifier, then we will prove the statement of the lemma. 

Step 1: From [38l Appendix B], |39l Lemmal] we get that there is a C > such that ||M^||L2(Tr2) < 
C||n„||2,2(T2). Additionally, Jg is a linear operator. [381 [39] also give \\Jef — /||l2(t2) < £||/II//i(t2) for 
/ G H^(T'^). Un is not regular enough to use this estimate, so we need a discrete version of this. 

As in the references above, using that 1 — < for E C, we find that 



(1-e- 



-£2|fc|2N2 



1 + 



< 



£2^2 if < 
5-2 if \k\ > 5. 



Hence, 



\uf, Mr), 11/^2(^^2-) — ^ ^ (l 6 



< ( sup 



(l-e-^'l^l')2 
1 + |A:|2 



<,2^ (l + |A;|2)|^(fc)| 



J2 (l + I^P) I'i^nik)]' 
2 



By Plancherel's/Parseval's identity we get immediately J2kez'^ 



\l2(^f2y Furthermore, 



l|-^^Afn'"n|li2(']r2) 



dx 



where Un (x) := Un{x + NnCj) for standard basis vectors ej, j G {1,2}. It's easily computed that 
ut'' {k) = S^(A;)e2''*^j/^", hence 

\\DN^Un\\l2^^2) = |S;(A;)|2[ Ar„(e2-'=i/^" - 1) ' + A^„(e2-'=^/^" - 1) 



Recognizing the difference quotient 

A;„(e2"^i/^" - 1) =27riA;i 

we deduce 



g27rifci/iV„ _ gO 

2mki/Nn 



2mki + 0{N-^), 



\DN,,Un\\\2(j2) =47r2 ^ |/c|2|m;(A:)|2 + C„ ^ |n;(A:)|2, 

k&I? k&I? 
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where C„ = 0{N^ ^). On the last term we can use again Parseval's formula. By the uniform bounds 
on 1 1 1 1^2(1^2) and ||-DAr^M„||x,2ci[-2) we then find that 

Ili2('ir2) < Ce'^, uniformly in n for n large enough. (B.l) 
Next we compute an estimate on the derivatives of uf^. 

and hence, since |nn(^)| < ll'iin||Li(T2)) 

d 



JeUn{x) < 27r||t;||ii(Tr2) ^ \ki\e " I'"' 



We compute 



and 



/"OO 

^ e-^''^' < 2 ^ e-^''^^ < 2 y e'^'^^ dfes + 1 = V^e^^ + 1 



fcaeZ fe2=0 



°° POO PC 



^e--'dx = e-^. 



fciez fci=o 



Because ||tin||Li(T2) ^ ll^^n||L2(Tr2) is uniformly bounded, we conclude (for e small enough) there is a 
C > such that \\VJsUn{x) ■ ek\\L°°{T'2) < Ce~^, € {1,2}. 

Step 2: Let 77 > and let n be large enough for the bounds proved in Step 1 to hold. Fix e > 
small enough such that, by (B.l), for each n we have \\un — '«nllL2(T2) < ^/3. 

By the bounds from Step 1 both ||ii^||^2(']r2) and \\VJ£Un{x) • CkWi'^iji-^ are uniformly (in n, for 
fixed e) bounded, and hence by the Rellich-Kondrachov compactness theorem |3H §5.7 Theorem 1], [U 
Theorem 6.2] the sequence {u^}^^^ converges strongly in L^(T^) as n — )• cxd. In particular it is a Cauchy 
sequence in L^(T^), so choose > such that for all n,m > Mg we have — "U^ 11^2(^2) < ??/3. 
Then for such n, m 

\\Un - Mm||L2(x2) < ||n„ - U^||l2(T2) + \\Um - UmWl'^iJ^) + \\Un " '"mllL2{T2) < V- 

Thus, {un}^=i is a Cauchy sequence in L^(T^) and therefore converges strongly in L^(T^). By the 
uniqueness of the limit it converges to u. □ 

Proof of Lemma 16.41 In both cases we assume without loss of generality that for is large enough 
such that ^ < I- For fixed iV > and i,j,k,l e 1^2, let x = {i/N,j/N) and y = {k/N,j/N). For 
z G T2 write /^'^'(z) = (/(x - z) - f{y - then from (O) 



(4,Jv).,>,/= E r''ir/N,s/N). 
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Also define Sl^n '■= {z € T"^ ■ \zi\ + \z2\ < L/N}. By repeated use of the trapezoidal rule for 
approximating integrals we then find 



L fL/N Q2jx,y 



Sl,n J~L/N CIZ2 



21 



s=—L ^ 

: / r^y{z)dz + RL,N, (B.2) 

'J St. AT 



'Sl 

where (Ci,C2) S 5'l,7V- By smoothness of / and compactness of we have \Rl,n\ < ^'fw^ some 
constant C/ > depending on /. 

For the first statement in the lemma we now find 

(47v)m,M =N^ I r'Hz)dz+N^RL,N = T^ [ r'^z)dz + N^RL,N 
JSl,n \^L,N\JSl,n 

4LV''^(0) = ^L^ifix) - f{y)f uniformly as iV ^ 00. 

The uniformity of the convergence follows by the bound of the smooth / on the compact domain T^. 
This proves the claim (since the composition of a continuous function and a uniformly converging 
sequence of functions is uniformly converging to the composition of the continuous function and the 
limit of the sequence). 

For the second statement the bound on / allows us to conclude that 

/ /""'^(z) / f'''y{z)dz uniformly as TV 00. 



The claim then follows by taking the limit — )• 00 in (B.2). □ 
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